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Abstract: The origin of the observed dark energy could be explained entirely within 
the standard model, with no new fields required. We show how the low-energy sector 
of the chiral QCD Lagrangian, once embedded in a non-trivial spacetime, gives rise to 
a cosmological vacuum energy density which can be presented entirely in terms of QCD 
parameters and the Hubble constant H as pA — H ■ mq{qq) /rrij^i ~ (4.3 • 10~'^eV)^. In 
this work we focus on the dynamics of the ghost fields that are essential ingredients of 
the aforementioned Lagrangian. In particular, we argue that the Veneziano ghost, being 
unphysical in the usual Minkowski QFT, exhibits important physical effects if the universe 
is expanding. Such effects are naturally very small as they are proportional to the rate 
of expansion iJ/Aqco ~ 10~^^. The co-existence of these two drastically difi^erent scales 
(^QCD ~ 100 MeV and H ~ 10"'^^ eV) is a direct consequence of the auxiliary conditions 
on the physical Hilbert space that are necessary to keep the theory unitary. The exact 
cancellation taking place in Minkowski space due to this auxiliary condition is slightly 
violated when the system is upgraded to an expanding background. Nevertheless, this 
"tiny" effect would in fact the driving force accelerating the universe today. We also derive 
the time-dependent equation of state w{t) for the dark energy component which tracks the 
dynamics of the Veneziano ghost in a FLRW universe. 
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1. Prelude 

It has been suggested recently [1] that the solution of cosmological vacuum energy puzzle 
may not require any new field beyond the standard model. The idea is based on the phi- 
losophy that gravitation can not be a truly fundamental interaction, but rather it must 
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be considered as a low-energy effective quantum field theory (QFT) [2]. In such a case, 
the corresponding gravitons should be treated as quasiparticles which do not feel all the 
microscopic degrees of freedom, but rather are sensitive to the "relevant excitations" only. 
In this framework it is quite natural to define the "renormalised cosmological constant" to 
be zero in Minkowski vacuum wherein the Einstein equations are automatically satisfied as 
the Ricci tensor identically vanishes. Thus, the energy-momentum tensor in combination 
with this "bare cosmological constant" must also vanish at this specific "point of normal- 
isation" to satisfy the Einstein equations. From this definition it is quite obvious that 
the "renormalised energy density" must be proportional to the deviation from Minkowski 
spacetime geometry. With this definition the effective QFT of gravity has a predictive 
power. In particular, it predicts the behaviour of the system in any non-trivial geometry 
of the spacetime. 

The first application of this proposal was the computation of the cosmological constant 
in a spacetime with non-trivial topological structure [1]. It was shown that the cosmological 
constant does not vanish if our universe can be represented by a large but finite manifold 
with typical size L ~ H~^, where H is the Hubble parameter. The cosmological vacuum 
energy density pA in this framework is expressed in terms of QCD parameters for Nf = 2 
light flavours as follows: 



This estimate should be compared with the observational value /ja ~ (2.3-10~^eV)^ [3, 4, 5]. 
The deviation of the cosmological constant from zero is entirely due to the large but finite 
size L of the manifold. 

The non- vanishing result (1.1) can be understood as a Casimir-type of vacuum energy 
when the boundary conditions and topology play a crucial role. The estimate (1.1) is 
based on our understanding of the ghost's dynamics: it can be analytically computed in 
the 2d Schwinger model and hopefully it can be tested in 4d lattice QCD as explained 
in [1, 6]. Moreover, one can fix the value of L today by equating the observed pA to (1.1). 
This prediction could be directly tested in future cosmic microwave background (CMB) 
data [7]. 

This contribution to the vacuum energy is computed using QFT techniques in a non- 
expanding universe without gravity. As it stands, it cannot be used for studying the 
evolution of this vacuum energy density with the expansion of the universe. In order to 
do so one needs to know the dynamics of the ghost field coupled to gravity on a finite 
manifold. In this paper we deal with the dynamics of the ghost field when eq. (1.1) is 
recognised as the strength of its potential energy. We will derive the dynamical equations 
which govern the evolution of the Veneziano ghost [8] in a curved spacetime. The analysis 
of these equations gives us a time-dependent equation of state w{t) which can be compared 
with observations. 

It is also interesting to note that a somewhat similar estimate was given in 1967 by 
Zeldovich [9] who argued that p\ ~ nip/Mp: this is numerically of the same order of 
magnitude as (1.1) if one replaces — )• H. Since then, the form HAqqj-^ has risen its 
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head several times, see for instance [10, 11, 12, 13, 14, 15]. Despite the apparent similarity 
between the form HAQ^jy and eq. (1.1), it is important to notice that, first of all, in our 
case it is the inverse linear size of the embedding compact manifold which actually appears 
in eq. (1.1) below, not the Hubble parameter, as in all the aforementioned papers. This is 
because our computations are done within QFT on a compact manifold size L, while we 
use H ~ in our numerical estimates. Therefore, our result is of fundamentally different 
origin. Finally, let us stress that equation (1.1) is not an estimate based on (although very 
well motivated) physical arguments, but a result of a precise calculation, which can be 
explicitly worked out completely analytically in a simplified 2d model defined on a finite 
manifold such 3jS 3i torus [6]. 

The paper is organised as follows. In the next section (section 2) we arrive at the equa- 
tions of motion which dictate the dynamics of the ghost field. In section 3 we discuss the 
interpretation of the obtained results using the standard lore of QFT which complements 
the analysis based on classical equations presented in section 4. In particular we argue that 
the physical Hilbert subspace changes when QCD is coupled to gravity. More specifically, 
we discuss how the Veneziano ghost being unphysical in the QFT formulation in Minkowski 
space exhibits non-trivial physics in a non-static universe. The corresponding effects are 
proportional to the rate of expansion H such that numerically they are naturally very 
small, iJ/AqcD ~ 10~^^. However, we claim that this "small" effect is the driving force 
which accelerates the universe now. In other words, acceleration from a vacuum energy 
with time-dependent equation of state w{t) is a direct consequence of the dynamics of the 
Veneziano ghost. Then, in section 4 we explicitly obtain the equation of state w{t) by 
solving the classical equations of motion for our fields. The classical solutions are analysed 
(mostly by numerical means) in section 5. In the same section we discuss observations 
and how they relate to our predictions. Section 6 provides some intuition about the scales 
involved in the mechanism, and their apparent "coincidence" (related to the eponymous 
problem in cosmology). We summarise our arguments in the concluding section 7. 

2. The ghost in Minkowski space: Lagrangian, gauge fixing, quantisation 

The mechanism that bridges the well known Veneziano ghost in QCD to the non-vanishing 
vacuum energy p\ ^ was recently put forward in [1], and aims at explaining the late 
accelerating phase of the universe as a result of a vacuum energy term which, in this 
scheme, is defined as the mismatch between fiat and infinite Minkowski spacetime and a 
topologically non-trivial one. 

If the cosmological constant is indeed a result of the spacetime we live in having a non- 
trivial topology, then there must be a messenger capable of carrying the information about 
the boundaries (and the related boundary conditions imposed on the quantum fields), and 
this is possible only when there are strictly massless degrees of freedom which can prop- 
agate at very large distances 0{H~^). We should emphasise that these are not necessary 
asymptotic physical states (in fact, they are not). Rather, they should be treated as a con- 
venient way to discuss large distance dynamics, non-trivial boundary conditions, invariance 
under large gauge transformations etc. Degrees of freedom of this kind are a handy way to 
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track the physics of the different constraints that arise when the introduction of a auxihary 
fields remove some non-locahties (anomahes) which always accompany gauge theories. 

The crucial for the present work observation is that while naively all QCD degrees 
of freedom can reach only very short lengths O(Aq^j-)), there is a unique (unphysical) 
degree of freedom which is exactly massless and can instead propagate to arbitrarily large 
distances: this is the aforementioned Veneziano ghost [8]. In short, the Veneziano ghost 
(see also [16]) which solves the U{1)a problem in QCD is also responsible for a slight 
difference in energy density between a finite manifold of size L and Minkowski space, 
such that p\ ~ HAqqy) ~ (10~^el^)^. Notice that this field is not in the physical spectrum, 
and as such it does not give rise to any of the usual problems associated with negative sign 
kinetic terms, propagators, norms and commutation relations. 

Thus, the correction to the vacuum energy density due to a very large (but not infinite) 
manifold is small, and goes as 1/L. The central point is that although small, this is 
not exponentially suppressed as exp(— L), as one could anticipate for any QFT where all 
physical degrees of freedom are massive (such as in QCD). Hence, the QCD ghost acts 
as a source for the cosmological constant, pA- This very small number (LAqcd)^^ — 
iZ/AqcD ~ 10~^^ nevertheless provides a non- vanishing cosmological constant which is 
surprisingly close to the observed value as eq. (1.1) explicitly shows. 

In this section we will write down the Lagrangian which contains all the dynamical 
information about the ghost field, and discuss how this field can be integrated out (as 
was done in the original paper [17]) to generate the rj' mass, which was the main result 
of [17]. If we were interested in rj' physics in Minkowski space, this would be the end of 
the story. However, our target is different: we want to obtain these results in Minkowski 
space keeping an explicit ghost field in the system in order to generalise the construction 
to a curved spacetime. It turns out that the Lagrangian put together this way includes 
another paired massless companion along with the ghost field such that in 4d it identically 
coincides with the corresponding expression for the 2d Schwinger model [18]. Unitarity in 
this new description is not violated (and so are not all the usual founding properties of a 
sensible QFT) thanks to the appearance of the ghost's companion which exactly cancels 
all unphysical contributions in the same way as it does in the 2d Schwinger model [18]. 

Let us remark here that there will be no new physical results when the well known 
resolution of the U{1)a problem is simply reformulated in a different way. However, such 
a new treatment will be essential and incredibly handy in putting the system into a curved 
spacetime. In this case, as we shall see in the next sections, there will be a novel contribution 
to the vacuum energy arising from the ghost and its partner, as the exact cancellation 
between these two does not hold anymore. The entity of this failed cancellation will be 
proportional to the departure from Minkowski spacetime, i.e. of order H. 

2.1 The Lagrangian 

The starting point for our analysis will be the Lagrangian which includes the ghost pole 
as proposed first by Di Vecchia and Veneziano in [17] (a variation on this Lagrangian was 
already available in [19]-we thank J. Schechter for drawing this important reference to our 
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attention), whose general form reads 



£ = A + Id.v'd^v' + -(0-^]q (2.1) 



+ Nfmq\{qq) \ cos 



rj' 

fv' 



where we exphcitly keep only the degrees of freedom we are going to be working with, such 
as rj' and the topological density q, while all others (including vr, K, rf) are assumed to be in 
Cq, and shall not be mentioned in this paper. In the rest of this section, unless explicitly 
stated, all quantities are defined in Minkowski spacetime. In this Lagrangian g.f. means 
gauge fixing term for three- form A^^p, see below, and the coefficient b ~ m,^, is fixed by 
the Witten-Veneziano relation for the topological susceptibility in pure gluodynamics. The 
topological density is defined as usual, 

q = ^ WG^^'^G"''^ = \^^^r.paG^'""' > (2-2) 

with 

G jiupcr = dfj^Aypu dcrApiyp + dpA(ypi, di,Ap(jp , (2.3) 



A =^ 



AldpAl-A';,duAl-AldaA';, + 2gCa,,AlA'Al . (2.4) 



The fields A'^^ are the usual N"^ — 1 gauge potentials for chiral QCD and Cabc the SU{Nc) 
structure constants. The constant f^' — f-n is the r]' decay constant (/,r being that of the 
pion), while ruq is the quark mass, and {qq) is the chiral condensate. 

The three-form A^pu is an abelian totally antisymmetric gauge field which, under 
colour gauge transformations with parameter A" 



5AI = dpA'' + igCabcA'A';, , (2.5) 



behaves as 



Aj/pa ^ Ajjpfy -\- di/Kpcr 9pA[/(j da^pu 1 (2.6) 
hp, oc A^d^K- - Aldpk'' . (2.7) 

In this way the four-form Gp^pa is gauge invariant. The term proportional to 9 is the usual 
^-term of QCD and appears in conjunction with the r/' field in the correct combination as 
dictated by the Ward Identities (WI). The constant 6 is a positive constant which would 
give the wrong sign for the mass term of the scalar field q, the property which motivated 
the term "Veneziano ghost" . However, this positive sign for h is what is required to extract 
the physical mass for the rj' meson, m^, ~ 6, see the original reference [17] for a thorough 
discussion. 

One should emphasise that the gauge fixing term in (2.1) has to be not confused with 
the standard gauge fixing term for the conventional gluon field A^, as it is related to the 
fixing of the gauge for the three-form Apyp describing the Veneziano ghost and carrying no 
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colour index. One can interpret the field ^^^p as a collective mode which is represented 
by a specific combination of the original gluon fields, which in the infrared leads to a pole 
in the unphysical subspace. We know about the existence of this very special degree of 
freedom and its properties from the resolution of the famous U{1)a problem: integrating 
out the q field (as shown below) provides the mass for the rj' meson. 

One more remark about the coefficient b which enters (2.1), and which is the principal 
ingredient in solving the U{1)a problem. Its magnitude is determined by the topological 
susceptibility in pure gluodynamics (without quarks) as 

I d'x{T{q{x),q{0)})e=o = ^^b. (2.8) 

Of course b = to any order in perturbation theory because q{x) is a total divergence 
q = d^Rf^, see below. However, as we learnt from [8, 16], 6 7^ due to the non-perturbative 
infrared physics; in fact, rn^, ~ b. Since we believe this approach to be the correct one, we 
shall adopt it in our work. 



2.2 Integrating out the topological field q ~ e^^upaG^'^''" 

If one is interested only in the spectrum of the meson particles arising from the low-energy 
effective Lagrangian (2.1), it is convenient to integrate out the scalar field q by means of 
its equations of motion, since there is no kinetic term associated to it. This is indeed the 
procedure followed by Di Vecchia and Veneziano in their original paper [17], the outcome 
of which, as follows from (2.1), is 

£ = £o + kV9V-^fe-|^V (2.9) 



2 ^ ' ' 47V, V 
+ Nfmq\{qq)\ cos 



rf_ 

f-n' 



where all the dependence on the three-form Aypa has disappeared. This formula explicitly 
shows that rj' receives a non- vanishing mass in the chiral limit, m^, ~ b/2Nc 7^ due to 
the non-zero magnitude of the coefficient b, which enters (2.1) and (2.8). This formula 
also reproduces the notorious Witten- Veneziano relation for the topological susceptibility 
in pure gluodynamics if one substitutes b = INcin^, into the expression (2.8) for the 
topological susceptibility. 

However, since we are aiming at uncovering the details of the dynamics of the ghost 
field, we want to reproduce all the physical results reviewed above without integrating out 
the topological three- form A^^p describing the Veneziano ghost. As we already mentioned 
earlier, in a curved spacetime the Veneziano ghost couples to the gravity field; therefore, 
in curved space, we will not be able to get rid of the ghost, and we will have to deal with 
it explicitly. In order to simplify this problem, in what follows we shall separate (project 
out) the relevant longitudinal degree of freedom in three-form Ap^p] this degree of freedom 
is the one that contributes to the topological susceptibility (2.8). 
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2.3 Gauge fixing: finding the ghost 

With the grand scheme just outhned in mind, we shall now choose the Lorenz-like gauge 

{d^A^^p) ~ {d^K, - d,K^) = , (2.10) 

in which we will carry out our manipulations. It is the same gauge which was discussed in 
the original paper [17]. In this formula we have defined as 

i^M^W^'""' 9 = 9^^^^. (2.11) 

We choose to work only with the longitudinal part of the ET^ field because only this lon- 
gitudinal part determines the topological density q = dp,K'^, and eventually leads to a 
non- vanishing contribution to the topological susceptibility (2.8). Therefore, we write the 
longitudinal part of K^^ as 

K^ = d^<^, (2.12) 

such that the expression for the topological density takes the form 

q = d^K^' = n^, (2.13) 

where $ is a new scalar field of mass dimension 2. We notice that the gauge condition (2.10) 
is automatically satisfied with our definition (2.12). Now our Lagrangian (2.1) can be 
expressed in terms of the <I> field as follows 



C = Co + ]^d^r]'d''ri' + Nfmq\{qq)\ cos 



fv' 



(2.14) 



where we plugged in the coefficient b — t- 2Ncm'^, as the Witten-Veneziano relation requires. 
If we integrate out the field we return to the expression (2.9) which describes the 
physical massive r]' field alone. 

As usual, the presence of 4-th order operator <I>nn<^ is a signal that the ghost is present 
in the system and may be quite dangerous. However, we know from the original form (2.1) 
that the system is unitary, well defined etc, in different words, it does not present any 
problem associated with the ghost. The lesson we are about to learn will be pivotal in 
promoting the system to a curved background, when we will not succeed in integrating out 
the ghost due to the coupling to gravity. 

It is convenient to define a new field (f)2 which is a combination of the the original r]' 
field and $ as 

V' ^ (<A2 + , (2.15) 
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which serves to complete the squares in (2.14) in such a way that one can ehminate the 
term j d'^xij'O^ = — J d^xd^^d^rj' . The Lagrangian now takes the form 



^ = i^df,(j)2d''(p2 + Nfmg\{qq)\ cos 



$ [m^,a + □□] ^ . 



fr]' 



+ 



72" 



(2.16) 



It is now straightforward to repeat the known steps in coping with the higher derivatives 
term in (2.16), namely, one recognises that this operator hides an extra degree of 
freedom, and can be explicitly reduced in terms of the associated propagator as 



A, 



A, 



\\m.n 



p^o [Af {nij^' ,x) - Af {p, x)] 



(2.17) 



which is the sum of a massive m^i scalar, and a massless ghost-like scalar. This means that 
the ^ field corresponds to two degrees of freedom which is almost an obvious statement if 
one formally writes the inverse operator as follows 



1 



-□ 



-□ 



(2.18) 



In analogy with the 2d Kogut and Susskind (KS) model [18], we will call the massive scalar 
field as (p while the massless ghost is 0i. The final Lagrangian, explicitly including all the 
relevant terms, thus becomes 



C 



1 



(2.19) 



-m^,(^2 + Nfmg\{qq)\ cos 



+ 



fv' 



where all fields have now canonical dimension one in four dimensions. 

We claim that the Lagrangian (2.19) is that part of QCD which describes long distance 
physics in our context. There are no new fields or coupling constants beyond the standard 
model: the Lagrangian (2.19) is part of it. Notice that (2.19) is exactly identical to that 
proposed by Kogut and Susskind in [18] for the 2d Schwinger model, and therefore, our 
4d system (2.19) is obviously unitary, and satisfies all other important properties of QFT 
similarly to the 2d Schwinger model. 

As we have shown, the Veneziano ghost in QCD is represented by the (pi field in (2.19) 
and it is always accompanied by its companion, the massless field (p2- This is the basic 
reason why the Veneziano ghost in QCD does not lead to any disastrous consequence in 
the system as we shall discuss below. 

2.4 Unitarity and the ghost 

As in the KS model, we can readily proceed to quantise this theory, and in order to do 
so we will impose canonical equal-time commutation relations for the fields <j), (pi, and (p2 
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which follow from (2.19) as 



x) , dt(j){t, y) 
[(^i(t,x) , dt(t)i{t,y)] 



> {x-y), 
i5^{x — y) 
>^{x-y) , 



(2.20) 



whence we evince that (pi is a massless ghost field, and its propagator will have an opposite 
sign in comparison with the conventional fields. 

The cosine interaction term includes vertices between the ghost and the other two 
scalar fields, but it can in fact be shown [18] that, once appropriate auxiliary (similar to 
Gupta-Bleuler [20, 21]) conditions on the physical Hilbert space are imposed, the unphysical 
degrees of freedom (pi and (p2 drop out of every gauge-invariant matrix element, leaving 
the theory well defined, i.e., unitary and without negative normed physical states, just as 
in the Lorentz invariant quantisation of electromagnetism. Specifically, this is achieved 
by demanding that the positive frequency part of the free massless combination {(p2 — (pi) 
annihilates the physical Hilbert space: 



With this additional requirement the quantum theory built on the Lagrangian (2.19) is well 
defined in any respect, and the physical sector of the theory exactly coincides with (2.9) 
which was obtained by a trivial integrating out procedure [17]. Yet, there is one place 
where the ghost has physical consequences, and that is the mass spectrum which, through 
the Witten-Veneziano mass formula [8, 16] relates the mass of the r]' to the topological 
susceptibility of the model. It is precisely the topological susceptibility which enjoys the 
uncancelled contribution of the ghost (pi with its companion (p2 because the topological 
density expressed in terms of see (2.13), which is a combination of the ghost (pi and the 
physical massive r]' field (p, but not of the companion (p2. It is this property that leads to 
a non-cancellation of the topological susceptibility, and eventually to the -q' mass. 

To wrap up this subsection in one statement: we have shown that the Di Vecchia- 
Veneziano Lagrangian (2.1) is nothing else than the KS Lagrangian promoted to four 
dimensions, eq. (2.19). We are not claiming that we have produced a new result in QCD. 
Rather we describe the well known resolution of the U{1)a problem in a different way, with 
the ghost explicitly and safely present in the system. As we shall see in the next sections, 
such a novel point of view is of formidable help when we upgrade to curved space. 

3. The ghost in curved space: physical interpretation 

In the standard QFT set up the Veneziano ghost is an unphysical degree of freedom, and 
represents a massless pole in the unphysical subspace of the entire Hilbert space. As is 
known, notwithstanding its being unphysical, the ghost still has observable consequences for 
the physical subspace; for instance, it gives a non-vanishing contribution to the topological 
susceptibility and therefore solves the famous U {1)a problem in QCD as Witten [16] and 
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Veneziano demonstrated [8] (see the review [22] for a generic introduction, and ref. [1] for 
a discussion of some specific aspects linked to this paper). 

We shall argue below that once QCD is coupled to gravity the "would be" unphysi- 
cal ghost may lead to another observable consequence (in addition to the one mentioned 
above). Namely, the ghost field may acquire a vacuum expectation value in an expanding 
universe, such as our FLRW universe. All effects related to this "ghost condensation" are 
proportional to the rate of expansion H so that numerically they are naturally very small, 
H/Aqcd ~ 10~^^. However, this "small" effect is precisely the driving force which accel- 
erates the present-time universe. The effect of this "ghost condensation" in many (formal) 
aspects is very similar to the well known problem of particle emission in an expanding 
universe, see below, although with very different interpretation of the associated outcomes. 

Finally, we will see how collective modes of the ghost which can be excited in the 
expanding universe carry very small typical momentum ~ A; ~ i7 whilst higher frequency 
modes are strongly suppressed. Precisely this feature of the spectrum allows us to identify 
the extra energy with observed dark energy because the typical wavelengths of the 
excitations are of the order of the Hubble parameter, Afc~l/A;~l/i?~10 Gyr, and they 
do not clump on scales smaller than this, in contrast with all other types of matter. In 
the limit for which H ^ 0, which corresponds to the usual Minkowski QFT formulation, 
all the modes move back into the unphysical subspace, and we return to the well studied 
case [8]. 

3.1 The Veneziano ghost in Minkowski space 

We take off from the relevant part of the Lagrangian for the Veneziano ghost and its partner 
as derived in section 2. This Lagrangian governs the dynamics of the massless ghost, (pi, a 
massless scalar filed (j)2 and a physical massive meson (p. 

1 ~ 1 1 "^1 ~ 



+ Nfmq\{qq)\J\f cos 



+ 



fv' 



(3.1) 



where J\f means normal ordering. 

Going to curved space one just needs to covariantise the derivative operators (which in 
fact, acting on scalars, are going to reduce to the usual partial derivatives). The partition 
function is given by 

Z = N j [V^] [V(j)2] [V^i] exp {iS} , (3.2) 
S = [ d^x^C , (3.3) 



where for Minkowski space one should put \/—g = 1. The bosons satisfy the commutation 
relations (2.20) from whose signs we learn that (pi is a massless (but interacting) ghost 
field, and its propagator will carry a negative sign (in Minkowski space) which is consistent 
with the negative sign for the kinetic term in (3.1). 
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As we have mentioned in section 2, the theory is well defined once we require that the 
positive frequency part of the free massless combination ((^2 — 4'i) annihilates the physical 
Hilbert space: 

{(1)2 - cbiY+^ \nphys) = . (3.4) 

This condition has been discussed in great details in the 2d case of the KS ghost [18], 
but it can be promptly generalised to our 4d Veneziano ghost as every step is a carbon- 
copy of the 2d model, only with different forms for the 4d and 2d Green's functions. Of 
course, this does not endangers any of the arguments presented in [18]. It is important 
to notice that the KS ghost in 2d and the Veneziano ghost in 4d enter the expression for 
the topological susceptibility alone, without their companions (j)2- The reason behind the 
topological susceptibility receiving a non- vanishing (negative) contribution from the ghost 
traces back to this property, and it is in agreement with the WI (in fact, the WI demand 
so). 

Let us remark here that neither Witten in [16] nor Veneziano in his [8] did discuss the 
ghost's dynamics because such complication can be entirely sidestepped if one considers 
only the gauge invariant sector of the theory. In fact, Di Vecchia and Veneziano in their 
original paper [17] integrated out the -F^jypo- field even without fixing the gauge; see section 2 
for details on how this procedure would be implemented in our set up. The resulting 
effective Lagrangian would be given by eq. (3.1) where now (pi = (j)2 = ^ such that only 
the massive (j) remains in the system. In this case this field corresponds to the usual t]'. As 
our goal here is not to study the mesons' mass spectrum of QCD, but rather to analyse 
the dynamics of (pi and 02, whose physics is much more involved in an expanding universe, 
once again we will instead (safely) neglect the dynamics of the massive physical (j) field by 
putting = 0. 

The subsidiary condition (3.4) which defines the physical subspace can recast as 

(flfc - bk) l^phys) = , (^phys|(4 - bl) = , (3.5) 
where we expanded (pi and 02 on a complete orthonormal basis Uk{t,x) and Vk{t,x) as 

01 (t, f ) = ^ \^akUk{t, x) + a\ul{t, x) , 

k 

(j)2{t,x) = '^^bkVk{t,x) + blvl,{t,x) . (3.6) 

k 

The equal-time commutation relations (2.20) are then equivalent to 

[bk, bk'] = , [bl bl] = , [bk, bl] = 5kk' , (3.7) 

for the 02 field, whereas for the ghost modes they satisfy 

[ak,ak'] = 0, [aj.,aj,,] = 0, [afe,a^/] = -4fc' , (3.8) 

where again the sign minus appears in these commutation relations. The ground state |0) 
is defined as usual 

afc|0) = 0, 6fc|0)=0, VA;. (3.9) 
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The sign minus in the commutators (3.8) is known to be carrier of disastrous conse- 
quences for the theory if (pi is not accompanied by another field (/)2 with properties that 
mirror and neutrahse it. As thoroughly explained in [18], the condition (3.4) or, what is 
the same, (3.5) are similar to the Gupta-Bleuler [20, 21] condition in QED which ensures 
that, defined in this way, the theory is self-consistent and unitarity (together with other 
important properties) is not violated due to the appearance of the ghost. 

To see this, one can check that the number operator N for and (p2 takes the form 

N = 5^ {bibk - 4afc) , (3.10) 
k 

while the Hamiltonian H reads 

R = Y,^k(blbk-alak) . (3.11) 

k 

With this form for the Hamiltonian it may seem that the term — a^a^ with sign minus 
implies instability as an arbitrary large number of the corresponding particles can carry 
an arbitrarily large amount of negative energy. However, it does not take more paper than 
a back of an envelope to check that the expectation value for any physical state in fact 
vanishes as a result of the subsidiary condition (3.5): 

(Hphys|H|?^phys) =0. (3.12) 

In different words, all these "dangerous" states which can produce arbitrary negative energy 
do not belong to the physical subspace defined by eq. (3.5). The same argument applies 
to the operator N with identical result 

(Hphys|N|?^phys) =0, (3.13) 

where we can see explicitly the pairing and cancelling mechanism at work. 

In the next subsection we shall see how this simple picture drastically changes when 
QCD is coupled to gravity. For our purposes it is sufficient to consider the system of the 
two fields (pi and (j)2 (which represent the Veneziano ghost and its partner, respectively) in a 
curved background, ignoring all other physical massive fields such as rj' , Nambu-Goldstones 
TT, K, all glueballs, etc. 

3.2 The Veneziano ghost in curved spacetime 

It is well known that there are inherent subtleties and obstacles when we attempt to 
formulate a QFT on a curved space [23]. In this case there is not a natural choice for 
the set of modes that on which the fields are expanded, these sets being closely related 
to a more or less "natural" coordinate system. Indeed, the Poincare group is no longer a 
symmetry of the spacetime and, in general, it would be not possible to separate positive 
frequency modes from negative frequency ones in the entire spacetime, in contrast with 
what happens in Minkowski space where the vector d/dt is a constant a Killing vector, 
orthogonal to the t = const hypersurface, and the eigenmodes (3.6) are eigenfunctions of 
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this Killing vector. The Minkowski separation is maintained throughout the whole space 
as a consequence of Poincare invariance. 

For our specific problem, i.e., the study of the ghost dynamics in a curved background, 
these considerations imply that there will be no simple formulation of the physical Hilbert 
subspace. Such a formulation in Minkowski space relies on the existence of a well defined 
positive frequency mode, which is in turn possible because the Killing vector d/dt is defined 
uniquely, since the Poincare group is a symmetry of the theory. A direct consequence of this 
is the fact that there are no privileged coordinates available in this set up and, therefore, 
there is no natural mode decomposition similar to (3.6). This means that a transition 
from a complete orthonormal set of modes to different one (the so-called Bogolubov's 
transformations) will always mix positive frequency modes (defined with the annihilation 
operators and bk) with negative frequency ones (associated with the creation operators 
aj, and foj,). As a result of this mixture, the vacuum state defined by a particular choice of 
the annihilation operators will not be "empty" once we switch back to the original basis 
defined by and bk- In other words, in curved space one should generally expect some 
relevant physical effects due to the ghost modes. In particular, as we shall see, they can 
give non-zero contribution to the energy (3.11). 

Such drastic, profound consequences arising in going from Minkowski to curved space 
should not be a surprise to anyone who is familiar with the problem of cosmological particle 
creation in a gravitational background, or the problem of photon emission by a neutral body 
which is accelerating. Only very few systems of this kind can be studied and solved exactly, 
see e.g. the reviews [23, 24]. However, the generic picture emerging from these analyses is 
amazingly simple: the transition from one coordinate system to another leads in general to 
non-vanishing Bogolubov's coefficients which mix positive and negative frequency modes. 
Eventually, it signals a physical production of particles stemming from the interaction with 
the gravitating background. 

The spectrum of the produced particles as well as the rate of production have been 
discussed in literature in great details. The most important outcome that is relevant for our 
work turns out to be the fact that the typical magnitude of the Bogolubov's coefficients 
is proportional to the rate at which the background is changing (the Hubble parameter 
H in case of an expanding universe, or the acceleration rate if we are studying photon 
emission by a neutral body), and to the total extent of this process, e.g., the total amount 
of expansion. 

The characteristic frequencies of the modes gravitaty can excite in this set up are of 
order of the Hubble parameter cok — H, whereas higher frequency ones are exponentially 
suppressed. This last result is easy to understand physically, because one expects the 
strength of the expansion to be able to excite modes for which ^ H, but not to possess 
enough energy to reach the higher end of the spectrum, that is, high k modes are only 
excited very inefficiently. 

Following the standard technique for the computation of particle production in a curved 
spacetime we consider, along with the expansion (3.6), a second complete set of-barred- 
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modes 



= ^ ^akUkit,x) +alul{t,x) , (3.14) 
k 

(j)2{t, ^) = ^ \bkVk{t, x) + blvKt, x) . 
k 

The new vacuum state is defined as 

afc|0) = 0, 6fc|0)=0, VA;. (3.15) 

Now, in order to study the new vacuum, we should expand the new modes and vj- 
in terms of the old ones. Following the notation of the textbook [23] we obtain 

«fc = ^ {akiui + I3kiu*i) , (3.16) 
I 

Vk = Y^ {a'kivi + f3'^,ivl) . 
I 

These matrices are called Bogolubov's coefficients, and they can be evaluated as 

an = {uk,ui) , f3ki = - iuk,u^) , (3.17) 
a'ki = ivk,vi) , P'ki = - ivk,vt) , 

where the brackets define the generalisation of the conventional scalar product for a curved 
space 

{^Pi,iJ2) = -iJ^Mx)'^f.r2V^dJ:^', (3.18) 

where dT,^ = n^dT, with a future-directed unit vector orthogonal to the spacelike 
hypersurface S and dT, is the volume element in S. Any complete set of modes which are 
orthonormal in the product (3.18) satisfies 

{uk,ui) = 6ki, {ul,ui) = -6ki , {uk,ui) = 0, (3.19) 
{vk,vi) = Ski, {vl,vi) = -5ku ivk,vf) = 0. 

Similar relations, of course, are also valid for the Uk and modes which appear in the 
alternative expansion (3.14). Equating the two expansions (3.6) and (3.14) and making 
use of the orthonormality of the modes (3.19), one obtains for the annihilation operators 



Ofe = ^ [aikai + Pikolj , (3.20) 
I 
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The Bogolubov's coefficients possess the set of properties 



{oilkoCnk - PlkPmk) " ^1 



(3.21) 



I 



{oilkPmk — PlkOimk) = , 



y^ f O^ifcCtmA: ~ f^'lkf^mk) — ^1 



I 



As one can immediately see from (3.20), the two Hilbert subspaces based on two possible 
choices of modes Uk and Vk, which appear in (3.6), and Uk and Vk, which instead enter 
in (3.14), are different as long as Pki 7^ 0,/?^^ 7^ 0. In particular, the expectation value of 
the Hamiltonian (3.11) of the k-th state in the barred vacuum (0|Hfc|0) is 



which is in sharp contrast with eq. (3.12), derived in Minkowski space. 

A remark is in order here. While a\.ak partakes in the expression for the Hamilto- 
nian with sign minus, it nevertheless gives a positive sign contribution to the expectation 
value as a result of an additional minus sign in the commutation relation for the ghost 
field (3.8). Hence, no cancellation between the ghost (pi and its partner (j)2 could occur in 
the expectation value (3.22), in net contrast with eq. (3.12). 

The best way to gain some insight on the result (3.22) is to assume that in the remote 
past and future the spacetime is Minkowskian. In fact, there is a simple 2d model with a 
specific profile for the expansion function a(t) interpolating between two Minkowski space- 
times which can be solved exactly. The outcome (see sec. 3.4 in [23]) is that, even in this 
plain example f3ki 7^ 0, which can be understood as a production of particles by the ex- 
panding background. In our case this should not be interpreted as actual emission of ghost 
modes, as the ghost modes are not asymptotic states in Minkowski spacetime, and there- 
fore they can not propagate to infinity; this constrasts with the conventional analysis [23]. 
More appropriately, one should interpret (3.22) as an additional time-dependent contribu- 
tion to the vacuum energy in time-dependent background in comparison with Minkowski 
spacetime. This extra energy is entirely ascribable to the presence of the unphysical de- 
grees of freedom. As we mentioned earlier, this is not the first time that a conventionally 
unphysical ghost contributes to a physically observable quantity. A clarifying simple 2d 
example on this matter can be found in [25]. 

In more complicated cases when in the remote past and future the spacetime is not 
Minkowskian, the formulation becomes very subtle as even a natural definition of particles 
is typically not available (the particle number will not be constant, a fact which makes 
its measurement inherently uncertain). We shall not discuss all these subtleties for the 
fields treated in the present work as they have exactly the same nature as for conventional 
particles, and they have been extensively discussed in the literature, see [23] and references 
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to the original works therein. Let us only point out that, while the number of produced 
particles may be a deceiving concept, and depends on whole prehistory of the spacetime 
for its interpretation to be sensible, the stress tensor is free from such kind of uncertainties, 
which is precisely the reason why we computed (0|Hfc|0), see eq. (3.22), rather than (0|Nfc|0); 
see also [25]. 

In order to understand the general pattern for the behaviour of the Bogolubov's co- 
efficients Pki in the expanding universe, it is useful to assume that the expansion is an 
adiabatically slow process. In this case the problem can be analysed using a WKB type 
of approximation [23] similar to the adiabatic approximation in quantum mechanics. In 
the context of quantum mechanics it is well known that if the rate of adiabatic expansion 
is very small (of order H), then all the transition matrix elements are proportional to the 
same rate H. Moreover, the transitions between states with very large difference in energy 
AE, are exponentially suppressed as exp{—AE/H), see e.g. [26]. 

In the case of QFT in a curved spacetime the technique which is required for such an 
analysis is much more involved than in quantum mechanics, but the basic reason for this 
suppression has the same nature. In fact one can argue that the means by which one studies 
the cosmological production of particles using the Bogolubov's coefficients as discussed 
above, are similar to instantaneous perturbation theory in quantum mechanics. Instead, 
the analysis of the same physical problem using adiabatic vacua is akin to employing 
adiabatic perturbation theory to calculate transition amplitudes in quantum mechanics. 
In any case, the aspect that is essential for the present presentation is the realisation 
that in all cases the high frequency modes will be strongly suppressed as exp{—uJk/H) in 
comparison with the expansion rate H. 

So, the basic result of this section can be formulated as follows: when QCD is coupled 
to gravity the "would be" unphysical ghost, although it still is not an asymptotic degree of 
freedom, nevertheless contributes to the vacuum energy density in an expanding universe. 
This time-dependent "ghost condensation" can, for the most part, be regarded as particle 
emission in an expanding universe, bearing in mind the subtleties reviewed in this section, 
i.e., no actual particles are being produced). All such effects are proportional to the rate 
of expansion H, hence, very small, as H/Aqqd ~ 10~^^, and are seen to be related to 
the tiny momentum ujk — k H available for efficient gravitational interactions, higher 
frequencies being exponentially suppressed. We have also seen that when we flatten the 
spacetime back to Minkowski, all these excited collective ghost's modes will be seized back 
in the unphysical Hilbert subspace. 

The chief ingredient in this discussion is the subsidiary condition (3.4) or (3.5), which 
ensures a peaceful coexistence of widely different scales: Aqcd ~ 100 MeV and H ~ 10~^^ 
eV. We shall elaborate further in section 6 on how a number of fine-tuning issues which 
most typically plague dark energy models [27, 28], possess simple explanations within the 
framework advocated in this paper. 

3.3 Ghost dynamics vs. Witten's approach 

The ghost we have been working with in this work, and whose activity is central our 
discussion, was postulated by Veneziano in the context of the U{1) problem. It is instructive 
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at this point to look at the different approach to that issue as formulated by Witten [16], 
where the ghost field does not ever enter the system. 

As long as we work in Minkowski spacetime the two constructions are perfectly equiv- 
alent as the subsidiary condition (3.4) or (3.5) ensures that the ghost degrees of freedom 
are decoupled from the physical Hilbert subspace, leaving both schemes with the identical 
physical spectrum. In a curved space, on the other hand, we argued that the "would be" 
unphysical ghost can produce a positive physical contribution to the energy-momentum 
tensor. The question arises naturally: where is the corresponding physics hidden in the 
language of Witten? 

To answer the question we start by recalling the crucial elements of [16]. The main 
object under focus is the topological susceptibility in pure gluodynamics defined as 

X=^|dx(0|r{Q(x),Q(0)}|0), Q^^G^G^'"^ (3.23) 

2 

We observe that x does not vanish in spite of the fact that the operator Q is the total 
derivative, and therefore x = in any order in perturbation theory. The next important 
element is the sign of x which is negative (i.e., it is opposite to what one should expect 
for a contribution from a conventional physical degree of freedom). This is precisely the 
reason for the introduction of the Veneziano ghost which would saturate the topological 
susceptibility with the required sign. 

On the other hand, Witten obtains the same contribution from the contact term which 
is represented by the equal-time commutator term [16] 

Xcontact = i Jd^x{0\ [K\x), Q(0)] |0) (3.24) 

= -(o|2g'»,(o)g^-(o)|o) 

(see also [29] for detailed computations). The point is that equation (3.23) along with the 
contact term (3.24) determine the dependence of the non-perturbative vacuum energy Eyac 
on 6 as follows from the definition x = ~ ^1°" 1 0=0 • Therefore, a non- vanishing value for the 
topological susceptibility x 7^ with a negative sign solves the U{1)a problem [8, 16, 17], 
(see also eq. (3) in [1] where the U{1)a problem is reviewed in this context). 

The next step is to consider the system in a curved background. It has been discussed 
in length in this paper how this generalisation can be performed using the ghost degrees of 
freedom. If we want to follow Witten's approach instead, the corresponding information 
does not go away, but rather it is hidden in equation (3.24) where now the correlation 
function must be considered in the more general background. Such correlation in a general 
curved space is not known, nor it is for our FLRW universe (it is the renormalised two-point 
function of a non-abelian, strongly interacting theory). 

In simple geometries such as de Sitter spacetime a few special examples are known, 
see [23] for an overview, [30] for the free scalar case, and [31, 32] for the abelian vector 
field case. In particular, a computation of the simplest expectation value of the scalar field 
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((/)(O)0(O)) at coinciding points in de Sitter background turns out to be a very complicated 
problem [30] which requires a delicate subtraction procedure on the Green's function in 
the given curved background. The resulting expectation value depends, of course, on the 
properties of the background field (gravity in this case) such as the rate of expansion H. 
These are the effects which have been treated in this paper using the Veneziano ghosts. 
Conversely, if we were dealing with the solution proposed by Witten, such physics would be 
shadowed by the properties of the Green's function (3.24) in a curved and time-dependent 
background and in the associated subtraction prescription. In Minkowski space this is 
simply a subtraction constant (as it was introduced by Witten [16]), which is just a number. 
Generalising to a curved time-dependent background this would be a "subtraction function" 
which, in principle, should be computable from (3.24) after a proper generalisation for the 
new spacetime. 

To conclude this subsection: we do not aim at solving this intriguing but prohibitive 
task in this work. Instead, we want to point out that the physical phenomena (described by 
the Veneziano ghost in our framework) do not disappear when we use a different approach 
as that advocated by Witten. The corresponding physics is hidden in the properties of 
equal-time commutators, which, we suspect, is much more difficult to study in a curved 
and time-dependent background. 

3.4 Detour to QED 

The main goal of this subsection is to demonstrate that a very similar effect related to the 
ghost degrees of freedom in a time-dependent background is not a specific feature of the 
Veneziano ghost in QCD, but is to be expected as a generic feature of any gauge theory 
when some unphysical degrees of freedom are present in Minkowski space, in particular in 
the well studied case of QED^. 

We start with the standard formulation of the Gupta-Bleuler constraint [20, 21] in 
conventional QED in Minkowski space: here the positive frequency part of the free massless 
combination d^j_A^~^^^ annihilates the physical Hilbert space 



As is known, this additional requirement ensures that {'Hphys\d^i-A^\T~Lphys) = 0. It is 
important to notice that the condition (3.25) is formulated exclusively for the positive 
frequency part; the negative frequency part is not mentioned in (3.25). In fact, one can 
not demand a stronger condition such as (S^A'^) J^phys) = as it can not be satisfied even 
for the ground state: the negative frequency half contains creation operators. 

When QED is quantised along these lines, all the problems related to the unphysi- 
cal polarisations e*-*^) and e^^^ are automatically straightened out in a way similar to our 
computations for the expectation value of the energy (3.12) ^. In different words, all 

^The conventional non-abelian Fadeev-Popov ghosts are much more comphcated objects and shall not 
be discussed here. 

^In the terminology of the present paper the ghost <j)i would correspond to the timelike polarisation 
e'°'-with negative kinetic term-while its partner 4>2 would be the longitudinal polarisation e*"^'. 




(3.25) 
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these "dangerous" states which can produce arbitrary negative energy do not belong to 
the physical subspace defined by eq. (3.25). All physical states are such that they contain 
a mixture of photons with polarisations e*-*^^ and e^^^ such that (3.25) is maintained in all 
gauge invariant scattering matrix elements. 

Let us now generalise this construction to a curved space (or even in the simpler case 
when the background is time-dependent). As is known, in this case there is no natural 
choice for the set of modes on which the fields are expanded. In general, it would be not 
possible to separate positive frequency modes dfj_A^~^^^^ and write a well defined prescription 
as in (3.25), at least not in the entire spacetime, just as in the Veneziano ghost example. 
This is a crucial difference in comparison with Minkowski spacetime when the vector d/dt 
is a Killing vector, orthogonal to the t = const hypersurface, and the eigenmodes in the 
expansion of the field A^^ are eigenfunctions of this Killing vector. The separation of 
positive and negative frequencies modes is maintained throughout the whole space as a 
consequence of Poincare invariance, and therefore, the Gupta-Bleuler condition (3.25) is a 
well defined constraint in Minkowski space in contrast with a time-dependent background 
when the condition (3.25) can not be maintained in the entire space. This means that a 
transition from a complete orthonormal set of modes to different one will generally mix 
different frequency modes, and as a result of this mixture, the vacuum state defined by a 
particular choice of the annihilation operators will not be "empty" once we switch back to 
the original basis. 

As for the Veneziano ghost, since a natural definition of particles in a curved back- 
ground is typically not available (the particle number will not be constant), see [23], the 
interpretation in terms of particles is ambiguous. On the other hand, the energy-momentum 
tensor is a well defined operator. It is our interpretation that in the stress tensor we obtain 
new contributions which are attributable to "would be" unphysical polarisations such as 
e^*^^ and e*^^-*, rather than claiming that these are actual (asymptotic) particles. Let us 
stress in passing that both e^"-* and e^'^-* contribute with positive sign to expectation value 
of the Hamiltonian, similar to eq. (3.22). 

As formulated above, this is not a statement about a possible violation of gauge in- 
variance. Instead, it is a statement about the global non-cancellation between unphysical 
polarisations in a time-dependent background. Moreover, we expect that in the "physical" 
Coulomb gauge Aq = 0, diA^ = (when one keeps only the physical degrees of freedom) the 
corresponding effects would also emerge. However, they will be hidden in the properties of 
the coinciding points Green's function in a curved background, similar to our discussions 
in the previous subsection 3.3. As is known, these correlators are quite non-trivial and 
very sensitive to the global properties of spacetime, in contrast with their simple behaviour 
in Minkowski space. We expect, however, that this new QED-related contribution to the 
dark energy would be negligible in comparison with what we have computed above for 
the Veneziano ghost, because QED effects, by dimensional reasons, are proportional to 
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4. Ghost classical dynamics and equation of state 



In order to do cosmology with our Lagrangian (2.19) we have to move the system to a more 
realistic universe which includes gravity. We want to analyse the dynamics in a curved 
background by studying the QCD Lagrangian which explicitly includes the Veneziano 
ghost (2.19), supplied by the auxiliary condition (2.21) selecting the physical subspace. 
As we demonstrated above, in Minkowski space these two elements faithfully reproduce 
the well known low-energy dynamics of (2.9). To this end we will first generalise the rele- 
vant piece of the QCD Lagrangian, given by (2.19), to a general curved background, and 
then try to solve the corresponding classical field equations for the expectation values of 
the fields. This is the standard procedure in studying the dynamics of scalar fields in the 
early universe, and we are making no exception here. 



4.1 Curved spacetime 

It is in fact merely a matter of rewriting the Minkowski Lagrangian (2.19) all what one 
needs to do to promote it to a curved space. Indeed, its simplicity allows the straightforward 
generalisation 



(4.1) 
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where the covariant derivative is defined as = 5^ + so that, for instance D^j^V^ = 

One can easily check that all the steps that led to the Minkowski Lagrangian (2.19) can 
be just as well worked out for a general background, see [6] for an example. Therefore, ig- 
noring the O-term contribution, the action for our system will be given by 5 = J d'^x^/—g £, 
where g = det g^u is the determinant of the metric. 

Accordingly, the energy- momentum tensor for the system of the (p, (pi, and (j)2 fields is 

' ' « ,4.2, 



V -9^9"'' 

= -D^cpDycf) - Dfj_(f>2D^(f>2 + D^cpiDyCpi + g^^C , 

which satisfies the energy conservation equation D^Tj^ = as it must. The equations of 
motion on which the stress tensor is conserved are cast as 

+ (f)2- <t> 
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Notice how the combination {(j)2 — (pi) is special in that it is a plane wave (in Minkowski 
space) : 

D^'D^{cl)2-(l>i) = 0. (4.6) 

We suggest that these equations of motion are the interesting part of QCD describing long 
distance physics in a curved spacetime. There are no new fields or coupling constants 
beyond the standard model. 

4.2 Solutions 

Let us now follow in some details the passages which lead to an (approximate) solution of 
these differential equations. From now on we will make the simplifying assumptions that 
the field (p sits at the bottom of its potential, and we fix its expectation value to zero. 
Moreover, we will be treating these fields only from a classical point of view (although 
we omit the angle brackets), and discuss the evolution of their vacuum expectation values 
(vevs) in the general background g^u , which for definiteness we choose to be of Friedmann- 
Lemaitre-Robertson- Walker (FLRW) type-see however appendix B for a generalisation 
involving compact 3d hypersurfaces (as required for consistency in our set up)-and hence 
given by 

ds2 = dt^ - a(t)2 ( dx^ + dy2 + dz^) , (4.7) 

where we denote the scale factor of the universe as a. For future reference, we define the 
Hubble parameter as H = a/a where the overdot is a derivative with respect to t. 

Let us change names for the two fields (pi (the ghost) and (p2 (its partner) to ^p and 
(p, respectively in order to emphasise that we treat ip V ^ classical fields in this 
subsection. Their quantum nature has been discussed in section 3. The classical treatment 
of the system is by far the widest chosen approach in dealing with dark energy matters, 
and this is the reason why we present this formalism here. However, while in other, 
perhaps more familiar, cosmological models such as infiation, the passage from quantum 
to classical is justified a posteriori (see for instance the discussion in [33], section 7.4.7, 
and reference to original works therein), and is a fundamental ingredient for the success 
of the inflationary theory, in coping with our quantum fields we do not expect such a 
"little miracle" to happen. In other words, the quantum nature of our fields which appears 
in their non-trivial dependence on the gravitational background as well as on the global 
properties of the manifold, is brought in at the level of the renormalisation procedure, and 
is therefore not describable in a purely classical approach. Therefore, all results based on 
the classical equation of motion presented in the next two sections should be taken very 
cautiously, as we do not have any theoretical arguments supporting the consistency of the 
classical approach (in contrast with the case of infiation). Still, we will make an attempt 
in this direction, these inherent difficulties notwithstanding, and address in more details 
the reasons for our mistrust in this naive construct thereafter. 

We want to work in conformal time r for the moment, which amounts to the coordinate 
transformation dr = At /a. Finally, we will initially solve for the rescaled fields and ip 
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defined as 



99 = a^p 



(4.8) 



Primes will mean derivatives with respect to conformal time 8/ dr. 

With these definitions one can either recast the action or act directly on the equations 
of motion. In any case, the explicit equations of motion for the tilded fields are 



where to simplify things we keep only the linear term in the expansion, s'm{(p2 — (pi) ~ 
{<p2 — (pi)- This approximation is more than sufficient for our qualitative discussions. The 
a-dependent mass-like term is defined as 



A few comments on this potential are in order here. First of all, the a" /a term is nothing 
else than the usual friction term in the original time t coordinatisation, for the original 
fields if and ip. This becomes an effective time-dependent mass term when one employs 
conformal coordinates and rescales the fields as in (4.8), and is proportional to the Ricci 
scalar curvature R. In our Lagrangian the fields have been chosen to be minimally coupled 
to gravity, but in general this coupling will appear and be of the form (S,^ip'^R — S,^ip'^R), 
and in the case of conformal coupling will make the a" /a term in (4.11) cancel. 

Second, the reader will have noticed that there is an extra factor of l/(m^/L) in the 
second term. This is a very subtle point which we have thoroughly discussed in our pre- 
vious paper [1], and it is a result of a subtraction procedure that compares the values of 
the vacuum energy in Minkowski space with that in a general compact manifold (such 
as a torus of size L). This prescription aims at extracting the physical and measurable 
portion of the vacuum potential energy of the ghost field, by taking the difference between 
the vacuum energy in compact space and that in infinite Minkowski space, and is of strict 
quantum mechanical origin. Essentially, this is our definition of the vacuum energy when 
the "renormalised energy density" is proportional to the departure from Minkowski space- 
time geometry and remains finite. This external (to the classical system) prescription will 
necessarily introduces some inconsistency in the classical picture, as we will see below. 

Notice that the correction (which was computed exploiting the topological suscepti- 
bility of QCD when the ghost is present) is linear in the inverse size of the manifold, not 
exponentially suppressed, as one would generally expect in the confined phase of QCD. 
This contribution is a direct result of the existence of the ghost and can be computed in an 
exact way in a compact 2d spacetime in the context of the Schwinger model [6] , and indeed 
shows the linear dependence we have used in (4.11). There are very strong arguments to 
believe that this is going to be the case in 4d as well, for the topological structure (with 
the contribution of the ghost to the topological susceptibility as demanded by the WI) is 




(4.10) 



(4.9) 
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identical in both models [1], and because as we have just shown in sec. 2, the low-energy 
Lagrangians coincide (2.19). There is therefore no arbitrariness in this choice of the po- 
tential: in fact, this is not a choice, as the potential is the one that correctly captures the 
chiral dynamics of low-energy QCD. 

Since ultimately we are interested in the combination {(p^ — ijP') which enters the 
expression for the equation of state, we look for solutions of the corresponding combination 
in conformal time for the rescaled fields. Even in curved space, the solution for {(p — ip) 
is straightforward. Assuming the fields do not depend on the three spatial coordinates we 
find 

^ ^(^-'0)=O (p-^f) = co + ciT , (4.12) 



with Co and ci the integration constants. The following step is to use this solution in the 
equation for (p (4.9) which becomes 

^" = -a;2 (co + cir) . (4.13) 

This equation can be solved analytically if we assume that R <^ Aqcd /L, as is certainly 
the case in the late-time regime we are interested in, in which case the curvature term in 
can be safely neglected. The time dependence of then simplifies greatly, being the 
same as that of a^/(a^L) oc a{t). To keep things simple we also choose ci = 0; of course 
the results we will present in the following section can only be generalised by retaining the 
full form for (p — ip). Equation (4.13) is then solvable for the first derivative of the field, 
and so is its counterpart for the original field p{t), whose solution reads 

Pit) = ^^-pH+ - '^co(t - to) , (4.14) 

where Iq is the present time, and obviously the conformal length L scales as L = LQ{a{t) / a{tQ)) . 
We can simplify this expression further, as long as we are working in the late-time regime 
(that is, much later the QCD phase transition epoch), by noticing that if the characteristic 
vev of the scalars 99 and ip is of order of the QCD scale (which we can parametrise as 
Aqcd or with /^/), although small enough to justify the simplification of the cosine inter- 
action term, then the first of all the terms in equation (4.14) will dominate. One implicit 
assumption we are also making is that the physical size of the manifold Lq is of order 
(but still larger) of the Hubble length, Hq, as required by observations, see the discussion 
in [7] (if this were not the case, then also the term w^(to)coto/a^ would give a similar 
contribution-notice that the dependence on time, through a, is the same). 

As mentioned earlier, we are aiming at finding the explicit expression for [p^ — ip'^), 
which leads us to the following expression 



p"^ - = (^(P + (p - V') 



Cl Co + ClT 

2p ^ H H 



Cl Co + cir 
~2 ^ 



-2p{tof-^, (4.15) 
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according to the simplifications explained in the previous paragraph. 

As a last step, let us recall what cq is. In equation (4.12) we have defined it to be a 
constant in conformal time for the rescaled fields (p and if). This means that once we switch 
back to the physical fields, cq will not be constant anymore, but is now proportional to 
a(t), and in our normalisation where the initial conditions are given today, it is going to 
be proportional to 00^3(^0) — Q^o/jj'- This is going to be the case for the other contribution 
coming from (/j(to) which is going to be of order /^,. All together, since we do not know 
the position of the field in its potential now, nor its velocity, it is simplest to parametrise 
them by means of one and the same constant (which we therefore expect to be somewhat 
smaller than 1) baptised c^p ~ ^{to)^{tQ) / f^, (in principle positive or negative). 

The final expression for the kinetic term for the two fields f and ip (or (f)2 and (pi, 
which is equivalent) is simplified to 



where we have defined ai = 2c^f^,Ho. A few remarks on the estimate (4.16) are in order. 



there are no new fields and coupling constants involved in our analysis, as everything is 
within QCD. The unknown coefficient c^p ^ 1 which appears in our description is not a 
consequence of unknown physics. Rather, it is determined by the initial conditions when 
ip{t) and ip{t) are treated classically. 

Lastly, the result (4.16) may look suspicious as it obviously violates the subsidiary 
condition (2.21). Let us remind that this condition was a key element leading to the 
cancelation of the ghost contribution (pi with its companion (p2 such that a QFT based 
on (2.19) is well defined in all respects as explained in section 2.4. While the condition (2.21) 
is indeed violated in a curved space, we should notice that the level of its violation is 
proportional to the rate of expansion H such that in a Minkowski spacetime, when H = 0, 
we again obey (2.21). Moreover, according to the logic of this paper, the breaking of 
the requirement (2.21) in curved space, which ultimately is proportional to Aqqy)Hq, is 
part of the energy density that accelerates our universe today. Why and how the physical 
subspace changes when the system is defined in a curved background has been explained 
in depth in section 3, where we explicitly demonstrate that the vacuum expectation value 
of the Hamiltonian representing the ghost field (and its companion) in such set up does 
not vanish. We shall argue that it is instead proportional to the expansion rate H, which 
is consistent with our estimate (4.16) based on the analysis of the classical equations of 
motion. 

4.3 Equation of state 

We have finally come up with an analytic and well defined expression for the sum of 
the kinetic energy densities of the two fields and ip, eq. (4.16). Their total energy 




(4.16) 



First, as expected, this contribution is the same order of magnitude, Aqq-^Hq, as the 
potential term, Aqqj-,/Lo, discussed previously [1]; both contributions are of the same order 
of magnitude of the critical density pc = SH^Mp. Second, as we underlined previously. 
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density (comprising the potential term as well) is found as usual by inspecting the energy- 
momentum tensor (4.2), which for a perfect homogeneous and isotropic fluid can be written 
in our conventions as 




{p + P) u^Uy + Pg^ 



(4.17) 



where p is the energy density and P the pressure density of some component; is the 
proper 4- velocity of the fluid. There are several subtleties related to the use of this perfect 
fluid form in the context of this work: they are unfolded partly in section 5, and more 
thoroughly in appendix B; we shall employ this simple form for now. 

The total (kinetic plus potential) energy density and pressure density of the vacuum 
are found to be just as usual 



with the second (still to be determined) constant appearing above being defined as 02 = 
(A'^^m^/^)/(4m^'Lo) (we have rewritten the chiral condensate in terms of measured pa- 
rameters TTT-TT and /tt as 'mq\{qq)\ = m^/^/4). 

The two expressions (4.18) and (4.19) can be combined in u) = Pa/pa to yield the 
effective equation of state w of the system with the Veneziano ghost and its partner (which 
is identified with the equation of state for dark energy), and which can be written as 



The equation of state (4.20) possesses a highly non-trivial time dependence as H in the 
r.h.s. of the equation implicitly depends on w{a). As for the well known scalar field this 
equation of state oscillates between -1 and 1 as long as the constants ai and 02 have the 
same sign (see however the discussion in section 5 and appendix B). We will be able to 
find solutions for a negative value of this ratio as well, which would imply an equation of 
state w < —1, which implies a phantom-type of field, and which could be useful in the 
construction of bouncing cosmologies (see the review [34] with an extensive references list 
of the original proposals). 

We are going to study this equation of state in section 5, explicitly solving it by 
numerical means within the redshift range < z < 2, which is the range spanned by 
direct measurements of the Hubble parameter, and therefore is the most interesting and 
constraining for the comparison of the model against the data. 

A final comment before moving on to the numerics of the Veneziano ghost. As we 
mentioned earlier, the numerical value of the constant ai which triggers its late-time mag- 
nitude is of the same order of magnitude as the potential part, see (5.5) below. This is why 
one should be careful in using the estimate (1.1) to directly infer the equation of state of 
this vacuum energy in an expanding universe, which would be mistakenly deduced to be 
w = —2/3 (corresponding to a~^), while it is actually determined by (4.20). Furthermore, 
there is one further subtlety in this respect, which is elaborated on in section 5. 




(4.18) 



(4.19) 



'w{a) 




(4.20) 
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5. Numerical analysis: confronting with observations 



5.1 Background 

The reference equation in this section is the definition of the time-dependent equation of 
state (4.20). This equation is highly non-trivial to solve, because the equation of state w{a) 
appears in the r.h.s. as well, being hidden in the explicit expression for H. Moreover, since 
in general w depends on a as we have seen, w = w{a), it appears inside the integral which 
gives the evolution of the vacuum energy density with time. 
The Hubble parameter is given by 



(5.1) 



The density fractions Q^j and il^ are today's values of JIa/ = Pm/Pc and = pk/ Pc, 
respectively. Here pc = 3MpH^ is the critical energy density for which the 3d space is 
flat. The suffixes M and A refer to matter (the sum of visible and invisible cold matter, 
mostly in the form of dark matter) and vacuum, respectively. The Hubble parameter as 
given in (5.1) comes from the Friedmann equation 



P 



2Ml 



K 



(5.2) 



with K the 3d curvature, and p is the sum of all the energy densities in the universe. Hence, 
3d flat space corresponds to k = which implies = 1 where, analogously, ^ = ilj (i 
indicises the different components of the universe). 

For the purpose of the numerical study we are going to present it is particularly conve- 
nient to normalise the scale factor to today's value oq, and we will also similarly normalise 
the Hubble parameter to Hq. In this case one can rewrite the equation of state (4.20) as 



w\a] 



1/2 



f)0^a-3 + f]0e-3/da(l+«.)/a 



1/2 



(5.3) 



aa'- 



where the constants a\ and and their ratio a are defined as 



Arrir 



-CgHo, a = 0:2/01. 



(5.4) 



In the preceding definition for 02 we follow [7] and compare the unknown size of the 
manifold Lq to the Hubble parameter today Hq via Lq = [cgHo)"^, with Cg an unknown 
constant of order one. Let us stress once more that this unknown constant, as was for dp, 
is not a consequence of unknown interactions of unknown fields. Rather, Cg is determined 
by the size of the manifold we live in as explained in [1, 7]. As one can see from (5.3) the 
equation of state w{a) is only sensitive to the ratio of these two unknowns and Cg. 
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5.2 Subtleties 



A few comments on equation (5.3) are in order. As the reader will have noticed, we are 
using an expression for the potential energy of the two fields (pi and (j)2 that comes from a 
very specific subtraction procedure defined in compact non-expanding spacetime [1]. There 
are two subtleties of this prescription which affect our modelling. 

a). First of all, the subtraction in Minkowski spacetime gives the potential energy an 
inverse linear dependence on the size of the manifold we started out with. This is fully 
analytically trackable in 2d [6], and is certainly present in non-expanding 4d QCD as we 
argued previously [1]. However, the system we are dealing with in this paper is embedded 
in a much more complicated set up, where gravity now partakes in the dynamics. It is licit 
to ask what would be the result of the same subtraction prescription when applied to a 
non-static case. 

Ideally, one would write down the most general form for the propagators of our fields 
, 4>2 and (f) in the general FLRW background in the compact space, with the most general 
dependence on the parameters that define the basic manifold (known as the Teichmiiller 
parameters, see appendix B). Once the propagators have been computed for compact 
and infinite FLRW spacetime, the subtraction procedure should be implemented, and the 
result would be the exact form of the "physical" potential energy density as defined in our 
framework. Notice that our "point of renormalisation" in the end will be Minkowski space, 
in such a way that the subtraction is performed in two steps, first in one and the same 
spacetime extracting the physical Casimir energy, and second renormalising to Minkowski 
space. 

Unfortunately this is a stunningly complicated task, and we are not able to tackle it 
analytically. However, it is easy to guess which kind of form the exact expression in a 
compact FLRW spacetime will be, or at least, what its principal contribution would be. 
Indeed, there are few dimensional parameters describing our system: the Hubble parameter, 
the sizes which define the manifold (of which we pick a characteristic size Lq), and QCD 
parameters of order /^/. In general, starting from the original expression 

l/ghost - ^ , (a-f J 

the corrections due to the expansion would be proportional to some dimensionless com- 
binations of the aforementioned parameters, which basically means that eq. (5.5) will be 
multiplied by some function of x = HL ~ 1 as the two other possible combinations H/ f^i, 
or {Lf^t)~^ are very small, O(10~^^), and can be safely neglected. Therefore, we phe- 
nomenologically introduce a function / with the intent to parametrise the outcome of the 
"ideal" subtraction: f{x) = f{HL). Remember that H and L are both normalised to 
present day values. 

The exact form of this function f{HL) is unknown, and should be fixed by the fully 
consistent calculation of the propagators for the ghost and its partner in compact FLRW 
universes as described above. In this work we shall limit ourselves to impose some physical 
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conditions on this function which would allow us to write down a representative form for 
it without actually knowing it. The requirements this function must meet are: it has to 
reproduce our Minkowski findings when the curvature is zero, that is, /(x) — >• 1 when 
X — >• 0; second, when the manifold size L — t- oo the dependence on the size itself in the 
energy density must disappear, which means that in the limit x — >• oo we ask that f{x) — t- x. 
A function which is compatible with these demands is 

f{x) = l + x, x = HL, (5.6) 

and the potential is accordingly modified in everywhere it appears (most importantly, 
in (5.3)) as 

^ghost = -^^f{HL) . (5.7) 

We should point out that the introduction of f{HL) is a phenomenological way to account 
for the unavoidable mismatch that is general would arise between the equation of state 
determined by our equations of motion in the FLRW background, and the simple promotion 
of the Minkowski result oc 1/L to the new background, having treated the potential by 
means of our subtraction procedure in a non-expanding universe only. 

b) . There is a second subtle point in the expression (5.3), namely the fact that it has 
been derived assuming a perfect isotropic fluid's stress tensor, while it has been pointed out 
by us previously [1, 7] that we need to work in a generally non-isotropic compact manifold. 
For instance the three-torus can easily be not isotropic, as its physical fundamental sizes 
can be all different, and the expansion rates be unequal in each direction, see appendix B. 
The perfect fluid form may not be justified for use in a sufficiently nontrivial manifold, and 
this will introduce more parameters in the equation of state, which for simplicity are not 
discussed here. Nevertheless, this introduces yet another source of error in our numerical 
solutions shown below. 

c) . The arguments outlined in the previous section about the validity of the classical 
analysis apply here, and inevitably lead to inconsistencies in some of the formulas we will 
be employing below. 

All this should be kept in mind in playing with the seemingly innocuous equations 
describing the system. The missing information about the manifold is going to appear in the 
numerical results in forms of inconsistencies between different parametrisations of the same 
quantities which ordinarily would be equivalent. As an example, the main discrepancy is 
going to be found in the expressions for the energy density, which can be directly calculated 
using the classical solutions of the equations of motion eq. (4.18) or by parametrisation 
through the effective equation of state found in (4.20). These subtleties are the topic of the 
appendices, and we postpone the details until then. Dealing with these uncertainties, and 
minimising their impact results in fixing the only two free parameters in this formulation, 
that is C(j, and c^: the precise way this is done is explained in appendix A, and here we 
only employ the outcome of that analysis, that is 

^ 1.1 • 10"^ and Cg w 5.1 • 10"^ . (5.8) 
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Recall that C(f, can be interpreted as the position of the ghost field in its potential times its 
velocity at the present time. This interpretation, and the numbers given above, are con- 
sistent with the approximations we have made in deriving the solutions of the equations of 
motion. Finally, Cg parametrises the size of the manifold in terms of the Hubble parameter 
today, which, for (5.8) means that L ~ l^H^^ , and consequently safely beyond the last 
scattering surface [7]. 

5.3 Results 

We plot in figure 1 the equation of state w{a) as a function of the scale factor of the 
universe, when the c constants are fixed as in (5.8), for the redshift range < z < 2, where 
the redshift is defined asl + z = l/o (remember that we normalise everything today). The 
numbers we adopt in our numerical analysis, beside (5.8), are: 

m^ = 135MeV, i^o = 2.13/i • lO^^'^GeV , 
m^/ = 958 MeV, /i = 0.71, 

= 133 MeV, = 0.73, ^ " ' 

fr,' = 104 MeV, Qm = 0.27. 

The equation of state describing the ghost's (and its partner's) dynamics depends on 
time as expected, and its today's value is 

w{z = 0) = Wo ^ -0.75 . (5.10) 



Equation of State 

2.2 1.8 1.4 1. 0.8 0.6 0.4 0.2 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Scale Factor / Inverse Redshift 

Figure 1: The equation of state solution of equation (5.3) in the text. The upper x-axis is the actual 
redshift. The equation of state is redshift (i.e., time) dependent, and today it takes the value wo « —0.75. 
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Two facts should be mentioned concerning this result. First of all, it is a representative 
result and while it is self-consistent, it should be not considered as the final word of this 
framework due to a number of subtleties discussed above and more deeply in appendix A. 
Had we been able to compute the propagators in the compact manifold and perform the 
subtraction as prescribed in an entirely analytical way, the will still be freedom of variation 
of these two parameters. In addition, the function /(x) was introduced to supply for 
our inability to compute the required Green's function on a non-static compact manifold. 
While its asymptotical behaviour is predictable, a precise expression for f{x) is not known. 
Our phenomenological parametrisation is not unique, and different fs conduce to similar, 
but still different results. It is not our intention here to play around with different choices 
of f{x) satisfying the same asymptotic behaviours; we instead want to present a typical, 
representative result with the simplest function that respects the imposed asymptotics. 

Second, the equation of state (5.10) is about 2a off according to WMAP [5]. However, 
this statement applies to an equation of state that does not evolve with time, and is 
significantly weakened when this possibility is taken into account. Indeed, WMAP in this 
case only measures the integrated equation of state from last scattering to now, which does 
not provide much information on its effective time dependence [35]. Moreover, the linear 
redshift dependence that is most widely used, that is, w{z) = wq + Wa(l — a) does not 
necessarily provide a good fit for the behaviour shown in fig. 1. 

One more comment concerns the sign of c^. Although it has been fixed greater than 
zero for it is the value that best meets the criterion enunciated before, we can also find 
solutions to eq. (5.3) in the opposite case. For instance, a solution with w{z = 0) ~ — 1 
approached from below (phantom) is available for \c^\ < 10~^ in the negative branch. 
This, as mentioned previously, could have interesting application in models of bouncing 
cosmologies [34]. 

Direct measurements of the Hubble parameter however can constrain the equation of 
state of dark energy directly, although such measurements are difficult to perform and in- 
herently bear large error bars. Using the most recent data available [36] we can graphically 
compare the standard cosmological model (ACDM, for Cold Dark Matter and cosmological 
constant) with the result of the equation of state we have just obtained. This is shown in 
fig. 2, where it is seen how the Hubble parameter does not fit the data points. 

The last graphical result we want to include is the behaviour of the different com- 
ponents of energy density of the universe for the same redshift range for which direct 
measurements of the Hubble parameter are available. In figure 3 we show the evolution of 
the vacuum energy component described in this work, together with the flat cosmological 
constant and the total matter energy density. 

Obviously, the vacuum energy density in our framework is no longer constant, and in 
principle it will meet that of ordinary and not-so-ordinary matter at earlier times than 
usually assumed based on ACDM model. However, the classical solutions we have found 
seem to point to an energy density which at early times grows faster than matter or 
radiation, thereby overwhelming them. With the choice endorsed in (5.8) we also find a 
model in which matter never dominates; a much smaller value for would have returned 
a better agreement with observations, although mining the consistency of the different 
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Predicted Hubble Parameter VS Data 
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Figure 2: Hubble parameter as a function of the scale factor of the redshift. The upper x-axis is the 
actual redshift. As explained in the legenda, the two lines correspond to the ACDM model (orange, or light 
grey in black/white) and the model proposed in this work (blue, or dark grey in white/black). The data 
are given with the corresponding error bars. 

expressions for the energy density (see below). For instance, we can obtain an equation 
of state virtually indistinguishable from a constant w ~ —1 in the desired redshift range 
foi" C0 ^ 10~^, which would re-establish a viable cosmology. A more serious problem 
appears to be the too fast growing of the kinetic energy density of the ghost in the past 
universe; this feature as well should not represent a too serious menace for this proposal, as 
it crucially depends on the assumptions we have made on the shape and time dependence 
of the potential. 

In general however, our results suggest that there will be appreciable amount of dark 
energy even at very high redshifts {z > 5) because of the specific features of the equation 
of state derived above. This is in contrast with the standard ACDM model where, in 
the same redshift region, the vacuum energy is completely ignorable. The study of this 
redshift span (z > 5) using Gamma Ray Bursts (GRB) as standard candles as suggested by 
ref. [37] (see also the analysis of [38] and references therein) represents a unique opportunity 
to verify or rule out our proposal. Indeed, as the mechanism described in this paper 
unambiguously predicts sizeable amounts of dark energy at very high redshifts, one can 
discriminate between models by studying z > 5 (see also the recent work [39] on constraints 
on dark energy contribution for z ~ (2 — 20). Finally, let us notice that a joint likelihood 
analysis of three different types of observations (supernovae type la data, CMB, and the 
Baryonic Acoustic Oscillations), similar to what has been lately performed in [40], is also 
sensitive to the peculiar features of our equation of state (4.20). 

Coming back to fig. 3, notice that there are two lines corresponding to the vacuum 
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Figure 3: Energy densities against scale factor. The upper x-axis is the actual redshift. As shown 
in the legenda, the solid blue (dark grey in black/white) hne tracks the vacuum energy density using 
expression (5.11), while the dashed hne follows (5.12). The flat black line is a constant cosmological 
constant, and the red (light grey in black/white) one is the energy density of total matter. 

energy because there are two possible ways of expressing it, either based on the parametri- 
sation through the effective equation of state 



with wla) as obtained graphically from the solution of equation (5.3) (see appendix A), or 
directly from the solution of the equations of motion for the ghost 



where f{x) is defined by eq. (5.6) and exhibits a very non-trivial dependence on z. As 
inferred from the plot, with the choice of c^j, and Cg in eq. (5.8) these two expressions 
are consistent within 15%. Let us stress that some mismatch between them (which are 
supposed to be exactly identical) is not a result of some fundamental flaw, but rather 
springs from our inability to compute f{x) analytically, combined with the fact that we 
are forcefully imposing a potential which has been derived in a purely quantum mechanical 
way (and that is a result of non-trivial boundary conditions) into the classical system. 
Therefore, these numerical results should not be treated as the final result: the full com- 
plete quantum calculation would return expressions that are automatically consistent one 
another, although their graphical representation by numerical means (recall that H and 
therefore w{a) itself appear in the expression for w = w{a) eq. (4.11)) will introduce once 
more some degree of numerical mismatch. 



PA(z) = PA(0)e-=^/^"('"+^)/' 



a 




Pa{z) = aiH{z){l + zf + a2(l + z)f{x) , 



(5.12) 
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5.4 The fate of the universe 



The evolution of the universe can be foUowed beyond today thanks to the fact that we 
have a precise equation tehing about the behaviour of the vacuum energy component at 
any time in the future. Indeed, one can look at the expressions for the equation of state, 
which, once a specific form for f{x) is selected, leads immediately to the evolution of the 
vacuum energy density. Looking deep into the future, the equation of state is going to be 



Thus, the equation of motion (for this f{x)) will slowly approach -1 from above, and the 
universe is dragged into a de Sitter state, just as in the usual ACDM model. 

Indeed, the energy density of vacuum will more and more approach a constant value, 
result that can be seen either in the parametrisation (5.11) when w{a) — t- —1 or in (5.12) 
since in this case the dominant piece will be the potential energy as given in (5.7) which 
when f{x) — )• x becomes proportional to H, which itself is approaching a constant. This 
shows that asymptotically the form chosen for f{x) is self-consistent. 

6. Fine tuning without fine tuning 

This section aims at explaining in a non-technical, intuitive way how a number of fine 
tuning issues such as "coincidence problems", "drastic separation of scales", "unnatural 
weakness of interactions", etc., which always plague dark energy models [27, 28], possess 
simple explanations within the framework endorsed in this paper, even without the need 
to introduce new fields, which come with new interactions, new coupling constants, and 
new symmetries. 

• We start with a close examination of the energy density scale pA ~ (2.3 • 10~^eV)^: 
where does it come from in our set up? Basically, this scale emerges as a result of our 
operational definition for the physical dark energy. Indeed, according to the scheme pro- 
posed in the introduction, the "renormalised cosmological constant" is set to be zero in 
Minkowski space wherein the Einstein equations are automatically satisfied. Thus, the 
energy-momentum tensor in combination with this "bare cosmological constant" must also 
vanish at this specific "point of normalisation" to satisfy the Einstein equations. From this 
definition it is quite obvious that the "renormalised energy density" must be proportional 
to the departure from Minkowski spacetime geometry. 

In the context of QCD on a compact manifold with characteristic size L, this extra 
piece gives precisely the appropriate scale for p\ as an offspring of known QCD dynamics 
at the QCD scale [1]. One should emphasise that this contribution is unique as all other 
vacuum contributions (such as gluon or Higgs condensates) vanish due to their cancellation 
with the vacuum energy defined in Minkowski spacetime, or are exponentially suppressed. 

• The second fine tuning problem goes under the name of "cosmic coincidence", or: 
why does it happen now? The same type of arguments which have been presented above 
suggests that the ghost vacuum energy should become relevant when its potential energy 




(5.13) 
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is of the same order of magnitude as pc- Equating these two quantities returns 

2Nf\ma{qQ)\ o o , . 

which can be read as an estimate for the age of the universe as to — -^o"^ — 10/i~^Gyr, 
which is indeed the correct hfetime for the present Hubble size. 

• Most of the proposals on how to attack the vacuum energy problem suggest to treat 
dark energy using a scalar field cp ("quintessence", "K-essence", etc., see the reviews [27, 
28]). The potential describing the dynamics of this new field (j) must be extremely fiat, 
with a typical mass jn^ ~ 10~^^eV [27, 28]. Such a tiny number is unheard of in particle 
physics, and certainly requires a great degree of fine tuning. 

In our model it so happened that the dynamics is also determined by interacting scalar 
fields, see section 2. However, while the interactions among them is strong (it comes from 
strongly interacting QCD), the potential energy which remains after the subtraction of 
the large Minkowski piece is indeed very shallow (6.1), and becomes even more shallow 
with time since the potential energy scales as 1/L ~ 1/a. The moral is that the effective 
potential is indeed flat, but this flatness does not stem from some approximate original 
symmetry (which is typically exploited in dark energy proposals [27, 28]), but rather is a 
consequence of the same governing principle stating that the physical energy density must 
be proportional to the deviation from Minkowski spacetime geometry. This difference in 
our case is 0{1/L). 

There are no new particles in our framework. Instead, there are the Veneziano ghost 
and its partner which are known low-energy fields of QCD and are responsible for the 
solution of the U{1)a problem. These fields are always present in the conventional QCD 
in Minkowski space. Their dynamics is trivial in Minkowski spacetime, and was there- 
fore ignored in all previous QCD applications. The scenario is analogous to that of the 
Gupta-Bleuler quantisation of QED [20, 21], where photons with time and longitudinal 
polarisations can be ignored in all S-matrix elements, because in the gauge invariant sector 
they always cancel each other. The ghost and its partner experience a new life in a cos- 
mological background, and lead to time-dependent vacuum energy in curved space QCD. 
They are still not physical (i.e., asymptotic) states. Nevertheless, their contribution to 
the vacuum energy is different from its Minkowski couterpart, and can be thought of as a 
time-dependent "ghost condensation". 

Finally note that QED photons, including unphysical polarisations, may also in prin- 
ciple contribute to dark energy. This contribution however is very small, as it is of order 
of or H'^ by dimensional reasons (see however the mechanism proposed in [41, 42]). 

• Not only well behaving scalar fields, but also ghost fields (similar to the Veneziano 
ghost with a negative sign for the kinetic energy term) have been extensively discussed in 
order to describe dark energy [27, 28] . Unfortunately, they normally lead to more questions 
than answers. As is known, ghost fields ("phantom fields" in the context of cosmology) 
are generally pestered by severe quantum instabilities, violating unitarity and other crucial 
fundamental principles of QFT. Why does our system not suffer from all these severe 
problems? The answer lies in the existence of the (p2 field, which always pairs up with 
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the Veneziano ghost (pi in QCD, see eq. (3.1). These two unphysical degrees of freedom 
and (/>2 drop out of every gauge invariant matrix element in Minkowski space, leaving 
the theory well defined, i.e., unitary and without negative normed physical states after the 
subsidiary condition (3.4) or (3.5) is imposed. 

As we argued in the previous section however, when QCD is defined on a curved 
background, the exact cancellation between (pi and (p2 does not hold anymore, and these 
fields do contribute to the energy (3.22). However, as we claimed, the departure from 
Minkowski space is controlled by 1/L ~ if rather than by the expected QCD scale Aqcd- 
This stems from the exact cancellation between (pi and (p2 that takes place in Minkowski 
spacetime. 

In particular, the constraint on the ghost maximal energy [kmax ^ 20 MeV) derived 
in [43] is obviously satisfied. Indeed, as we have already elaborated on previously, the energy 
which can be associated with the ghost in the expanding universe is only ~ A; ~ if. 
Therefore, since H ~ 10~^^eV <C kmax, we clearly meet the limit of [43] on possible 
dynamical violations of unitarity. Let us mention once more that this very low cutoff comes 
about because of the subsidiary condition (3.4) (or, equivalently, (3.5)) on the physical 
Hilbert space, and hence, it is a result of the dynamics of the ghost in curved space, rather 
than an ad hoc requirement on the form of the ghost's interactions. 

• Last, but not least, the same feature of the spectrum discussed above allows us 
to name this extra ghost dipole contribution (3.22) dark energy; indeed, since the typical 
wavelengths Xk associated with this energy density are of the order of the Hubble parameter, 
Afc~l/A;~l/i/~10 Gyr, the corresponding excitations are not going to produce any 
clustering on scales smaller than this, as opposed to visible or dark matter. 

7. Conclusion 

Dark energy has provided formidable headaches to physicists for over a decade. In trying 
to understand the physical principles and mechanisms, sapiently hidden by Nature, that 
are responsible for it, the physics community has produced a mole of models that is quickly 
climbing up to 4000 arXiv papers. On the experimental side, a great deal of efforts is being 
devoted in planning, projecting and performing tests which will help theorists pierce the 
veil of Maya of dark energy. 

In our previous paper [1] we have claimed that 4d QCD, once formulated on a com- 
pact manifold, will exhibit a very peculiar feature, that is, the vacuum state will depend 
linearly on the inverse characteristic size L of the manifold. This finding was supported 
by an explicit calculation in the context of the 2d Schwinger model [6], where the linear 
dependence can be tracked analytically all the way to the final result. The non-vanishing 
result (1.1) can be understood as a Casimir-type of vacuum energy when the boundary 
conditions and topology play a key role. 

The essential reason for this unexpected behaviour is to be traced back to the existence 
of a very special degree of freedom, the Kogut-Susskind ghost in 2d and the Veneziano ghost 
in 4d, which, despite being unphysical, are able to carry the information about the bound- 
aries of spacetime due to its protected masslessness. Once the existence and the form of the 
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QCD ghost's potential vacuum energy is established to be (1.1), the next logical step is to 
promote the background to a general curved background, in our case a FLRW spacetime, 
and study the dynamics of the system which is precisely the goal of the present work. The 
main achievements of these studies are listed below. 

1. The Veneziano ghost (pi is always accompanied by its partner (/>2. These two un- 
physical degrees of freedom (pi and (p2 drop out of every gauge invariant matrix element 
in Minkowski space, leaving the theory well defined, i.e., unitary and without negative 
normed physical states after the subsidiary condition (3.4) or (3.5) is imposed. The only 
place where the Veneziano ghost manifests itself is the topological susceptibility, because 
the ghost's partner (j)2 does not participate in its expression. However, it is through the 
topological susceptibility that we know about the existence of the ghost and its properties. 

2. This simple picture drastically changes when QCD is coupled to gravity, in which 
case another physically observable manifestation of the ghost emerges. In this case the 
"would be" unphysical ghost contributes a physical portion of the energy density of an 
expanding universe. All effects related to this contribution are proportional to the rate 
of expansion H such that numerically they are naturally very small, H/Aqcd ~ lO"'^^. 
However, this "tiny" effect is precisely what accelerates our universe today. In the limit 
for which — )• 0, which corresponds to the usual Minkowski QFT formulation, every ad- 
ditional contribution vanishes. We interpret this extra contribution as a time-dependent 
"ghost condensation" as explained at the end of section 3.2. We are not claiming that 
the ghost field becomes a propagating degree of freedom, or becomes an asymptotic state. 
Rather, we claim that our description in terms of the ghost is a convenient way to account 
for the physics hidden in the non-trivial boundary conditions as discussed in section 3.3. 
The physical phenomena (described by the Veneziano ghost in our framework) do not dis- 
appear when we use a different approach, for instance when the ghost is not even a part of 
the system, see section 3.3. However, we do not know presently how to track this physics in 
a ghost-free formulation. We strongly suspect that the corresponding description (if found) 
would be much more technically complicated in comparison with the approach presented 
in this work. 

3. The usual catalogue of fine tuning issues which typically pester us while tampering 
with dark energy does not appear here, as these questions possess simple explanations in 
our scenario, without the need to introduce new fields, which come with new interactions, 
new coupling constants, and new symmetries. The nature of this "fine tuning without 
fine tuning", see section 6, is not to be sought for in supersymmetry or any other extra 
symmetries imposed on the system (there are in fact, none), but it comes about from the 
auxiliary conditions on the physical Hilbert space (3.4) which accommodate the gigantic 
span of scales (Aqcd versus Hq). 

4. A word on the testability of our framework. There are in fact several places 
where one can look for supporting or ruling out evidence for our proposal. First of all. 
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as we have claimed in [1], the hnear dependence of the vacuum energy on the inverse size 
of the manifold can be tested by studying the topological susceptibility of QCD on the 
lattice by looking for its dependence on L. Second, we have argued in [7] that the size 
of the manifold is not overly large, and could in principle be tested in upcoming CMB 
maps (PLANK). Third, since in this paper we have provided the dynamical equations that 
govern the evolution of the ghost, its contribution to the Hubble parameter and to the 
expansion of the universe in general, its equation of state w{t) and its effects on structure 
formation, can all be confronted with cosmological observations (such as supernovae type 
la data, CMB, the Baryonic Acoustic Oscillations, Gamma Ray Bursts, see section 5.3 for 
details), a feature not often found in dark energy models. 

Finally, we would like to draw attention to the fact that it is quite amazing that QCD, 
which naively has no relation what so ever to the dark energy problem, may in fact be indeed 
the locus one should be scrutinising to uncover the dark energy secret, as we argued in this 
work. The emergence of the QCD scale in the vacuum energy may eventually unlock the 
mystery of the "cosmic coincidence" problem, that is, the observational similarity among 
visible and invisible scales: f^A — 417 a/ and Ojv/ ~ 5il.B where only 0,b (which represents 
the baryon contribution to fltot = 1) has an obvious relation to QCD, as the nucleon mass 
miy is proportional to Aqcdi while "God's particle" (the Higgs boson), contributes only 
few percents to ttiat (and correspondingly to Ob) through the quark's mass dependence 
on the Higgs vev. It is also fairly interesting that a number of dark matter puzzles (along 
with the dark energy mystery which is the subject of the present study), may also be 
related to QCD ^. The different sides of our construction allow for several types of tests to 
be performed, and some of them are radically new (and are entirely intertwined with the 
QCD nature of dark energy). Ultimately, observations will be able to pin down the winner, 
the model which most adherently follow the forms of nature. 
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A. Numerics 

This first appendix outlines some details of the numerical calculations which produced our 
figures. The first quantity we want to compute is the equation of state of the vacuum 

^We refer to few recent papers on the subject [51, 52] and references therein, where it has been argued that 
dark matter, in fact, may also have QCD origin. Tlie basic idea is that nuggets of very dense matter similar 
to Witten's strangelets [53] form at the same QCD phase transition as conventional baryons (neutrons and 
protons), providing a natural explanation for the similarity between Qm and Qb [54]. 



-37- 



energy fluid. This is given by equation (5.3) in the text which, once we incorporate the 
effects of the unknown function f{x) becomes 

where, following the main text, we choose f{x) = 1 + x and where x = HL. It is not 
possible to solve for ■w{a) analytically, as it appears both on the l.h.s. and on the r.h.r. 
under the integral 

f 1 + w(a) 

I = -3 da ^ ^ . (A.2) 



The easiest way to numerically overcome this problem is to assume that in the integral 
the equation of state is a constant, and then solve graphically for w{a). Of course the re- 
sulting equation of state will not be a constant, but it can be used back in the integral (A.2) 
which can be solved numerically and the resulting expression plugged back into (A.l) to 
obtain (once more, numerically) a better approximation for w{a). This is the procedure 
followed in obtaining the results shown in the figures in the text, figs. 1, 2, and 3. 

Practically, we have found that within two recursions we already have a reliable ap- 
proximation which does not change much employing further steps. The reason for this can 
be traced back to the fact that among the two components of energy and pressure density, 
that is, kinetic and potential term, at low redshift the potential term dominates (for our 
choice of the parameters) and the errors induced by mistaking w for constant in H mainly 
appear in the kinetic term. In the high z region, when the kinetic term dominates, there 
will be some more uncertainty, but this stabilises within a couple of steps. 

Once a plot for the equation of state has been obtained (which comes from feedbacking 
w into (A.l) once) we can use an interpolation of this to find the other quantities of interest, 
in particular the Hubble parameter and the energy density of the vacuum. 

As mentioned in the main text, since we are forced to introduce the function f{x) 
together with the fact that the information about the compact manifold is implicitly hidden 
in some of our expressions while it has been entirely discarded in others, will introduce some 
degree of inconsistency among different expressions for the same quantity. The clearest 
example of this can be immediately noticed by looking at the expression for the vacuum 
energy as given in equation (5.12): the potential term comes from a very specific subtraction 
procedure as explained in [1] and exhibits the 1/L behaviour which is one of the principal 
ingredients of our work. This subtraction has been defined and performed in a compact 
Minkowski space, and if promoted acritically to an expanding universe it would lead, in the 
case of a simple 3-torus with the same expansion rate in each direction, to the behaviour 
^ghost oc 1/a, which implies an equation of state w = —2/3, at least in the redshift range 
for which the potential energy dominates over the kinetic energy. At the same time the 
equation of state clearly shows that if this is the case then the equation of state approaches 
-1, in contradiction to what has been found just a few lines above. 

The resolution to this apparent paradox lies in the subtleties mentioned in the text 
and analysed here, that is, the expression (5.5) is based on our understanding of the ghost 
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dynamics without gravity: it can be analytically computed in the 2d Schwinger model and 
hopefully it can be tested in 4d QCD using the lattice QCD computations as explained 
in [1, 6]. This contribution to the vacuum energy is computed using QFT techniques in a 
non-expanding universe. As it stands, it can not be used for studying its evolution with the 
expansion of the universe. In order to do so one needs to know the dynamics of the ghost 
field coupled to gravity on a finite manifold, which we are not able to solve for analytically, 
and we therefore must resort to modelling it by means of the function f{x). 

Even doing so there will still be some inconsistency due to the fact that the information 
about the compact manifold that shows up in the potential energy indirectly influences 
also the resulting kinetic energy. However in this case the effects are subtler and it is 
not possible to easily model them as for the potential term. Our approach in this case 
is to use the freedom of choosing the two parameters c<^ and Cg to minimise the impact 
of these uncertainties on the quantities of interest, which will leave us with a consistent 
representative theory. This theory is not unique and may well not be able to fully reproduce 
the observational data, see the discussion in the main body of the paper. 

Thus, the operational criterion for the choice of the constants (5.8) is that the energy 
density of the vacuum computed using the effective equation of state or using directly the 
solutions for the equations of motion (formulas (5.11) and (5.12), respectively) must be 
the same within 15%. If one fixes the energy density of the vacuum today to = 0-73 
without errors, this automatically fixes one parameter in terms of the other as 

which in turn gives Cg in terms of c^. 

This criterion led to the choice (5.8), and will be different for a different choice of the 
unknown /(x). In figure 4 one can see how by moving away from the choice values the two 
expressions (5.11) and (5.12) progressively drift farer apart. Allowing for an error in the 
determination of the vacuum energy density today would open up some more parameter 
space where the criterion formulated above would be satisfied. Finally, notice that the 
resulting Cg is slightly different than the one inferred in our paper [7] but only marginally 
{cg ^ 0.059 in that paper, whereas Cg w 0.051 here), since at zero redshift most of the 
vacuum energy density is in fact stored in the potential energy, fact that confirms the 
validity of the hypothesis made (and the results thereby obtained) in that work. 



B. Compact universes 



The idea of a compact universe has surfaced in several places throughout the paper. It is 
only thanks to the embedding of the spatial hypersurfaces in a compact 3d manifold that 
we expect to observe a linear departure from the usual vacuum contributions in the form 
of or H'^. Moreover, although in a much more subtle way, the linear dependence of the 
potential energy as in (5.5) or (5.7) on the characteristic size of the compactified manifold 
also affects the solutions for the differential equations and therefore the kinetic piece (4.16) 
of the total energy density of the Veneziano ghost. 
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Figure 4: Difference in percent between the two parametrisations of the vacuum energy density for 
different values of c,^ as a function of the scale factor. As usual, the upper x-axis is the actual redshift. 



In the main body of the paper we have been rather generic in discussing this embedding 
in a compact manifold, the reason for this being the fact that in this way we could obtain 
interesting, self-consistent results with the minimal amount of free parameters (in fact, 
only one, c^). Consequently, taking into account the more general case of an anisotropic 
homogeneous expanding compact manifold would only enlarge the space of possibilities, 
however at the expense of clarity. Yet one more motivation for not carrying out the 
calculations in the most general set up is given by the acknowledgement that the most 
paramount quantity in our work, the value of the chiral condensate in a given background, 
is not known (and discouragingly hard to calculate), thereby eluding a completely closed 
analysis. 

For the sake of completeness, we will nevertheless review some of the key concepts and 
results in compact 3d geometries in this appendix. A clear and thorough treatment of the 
topic can be found in [44, 45, 46, 47]. Let us begin with a few definitions (in this appendix 
we employ the abstract index notation as in [48]). We are interested in "geometries" in 
the sense of a pair, {^^^ M , gab), where ^'^^M is a non-compact simply connected Lorentzian 
manifold, and gab its metric, or universal cover, a solution of the Einstein equations. This 
geometry needs to be reduced to a 3d Riemannian one that we call {M ,h{t)ab), at each 
given time (for the purpose of engineering a compact 3d slice of the 4d manifold the time 
lapse and shift functions can not depend on the spatial coordinates, but only on time itself), 
whose compact subgeometries can be studied. Here the term "subgeometry" refers to the 
fact that a geometry (X, hab) can always be thought of as a pair (X, H) where H is some 
symmetry group, and the corresponding i?-invariant metric can be always constructed: 
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therefore a pair (X,Q) with Q C H is a subgeometry of {X,H), and the same applies to 
the corresponding manifold- metric pairs. In this example H would be the group of the 
isometrics of the metric hab- 

In order to make the 3d manifold compact, and find the corresponding metric, we need 
to know the group of isometrics of {M,hab), which we denote IsomM (in fact, we need 
the subgroup of IsomM which is orientation-preserving and extendible) . Once this group 
is identified we make (M, hat) compact by taking its quotient with a discrete subgroup 
of IsomM, which we name T. This group is basically the group of the identifications 
that wrap and glue the manifold to make it compact, and it is referred to as the covering 
transformations group. Therefore, the resulting compact 3d geometry (M, hat) is metrically 
diffeomorphic to (M, hab)/^- At this point one can jump back to 4d and naturally extend 
the covering group T to the Lorentzian geometry (^'^^ M , gab)- the resulting compact 4d 
spacetime {^'^^ M, g^b) will inherit the metric from the universal cover gab, and if that is 
homogeneous, so will the compact space metric gab be. 

Notice that thinking of a compact space as of M /T automatically separates the degrees 
of freedom of M to those of T. The latter admits a group of smooth deformations (for 
instance, stretching or skewing a parallelepiped in 2d) which are known as the Teichmiiller 
deformations, which form a group and are parametrised by the Teichmiiller parameters. 
Therefore, the Teichmiiller parameters know about the global structure of the manifold 
while preserving all local quantities, whereas the local coordinate parameters are insensitive 
to the global structure. In general there will be degrees of freedom associated with the 
Teichmiiller parameters, and they may have some dynamical role, which for simplicity we 
neglect in our report. The complete classification of the possible compact topologies and 
their Teichmiiller parameters is given in [45, 46, 47]. 

With all the definitions at hand, we can present a specific example for the simplest 
compact manifold in a fiat (i.e., Euclidean) 3d space, that is, the three-torus T^. In 
general this type of compactification is admitted in the context of the Bianchi I or Bianchi 
VII(O) geometries only. These geometries of course admit also other more complicated 
orientable compactifications, e.g., T^/Z2 or T'^/X^ etc., each of which will come with its own 
Teichmiiller space. For concreteness, we stick to a non-twisted on Bianchi I geometry. 

The aim is to find the most general metric which contains all the information about 
the dynamical local variables and the global deformations of the manifold. Schematically, 
the process goes as follows. First of all, we identify the 4d universal cover, that is, we 
specify gab from which we immediately derive the 3d slice {E^, hab)- This is going to be a 
particular solution of the Bianchi I type, which for instance, once a coordinate basis has 
been chosen, is given by 

ds^ = dt^ - a{t)^ dx^ - 6(t)2 dy2 _ c{t)^ dz^ , (B.l) 

with three, in general unequal, scale factors. 

At this point one can choose a specific 3-torus Tg by giving its fundamental domain. 
The resulting geometry will therefore be {Tg,Hab) where Hab is the metric inherited from 
the universal cover hab (in fact, it is the same). However this is not yet what we want, 
as the degrees of freedom corresponding to the Teichmiiller deformations are hidden in 
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the definition of the compact manifold, that is, in its covering group. It is however always 
possible to transfer the information about the global structure of the manifold to the metric. 
This can be achieved by pulling back the pair {Tg,Hab) to (Tg , hat), where now Tq is the 
reference 3-torus (the unit cube), and hat contains all the global and local geometrical 
information. 

In order to calculate the metric hab we need the pullback from the torus Tg to the 
reference Tg one, which can be written explicitly in terms of the generators of the group of 
the isometries defining the covering group, which in this case are three dimensional rotations 
and translations, parametrised by the three three- vectors gi, g2, and (73. Precisely how this 
is done is detailed in [45, 46, 47], and we only repeat their results here. In the Bianchi I 
case given in (B.l) the result would be 



= a{tfglg] + }?{t)gh] 



citfgfg' 



(B.2) 



It is easy at this point to make the connection with the Teichmiiller parameters. Indeed, 
every metric on the reference torus Tg can always be parametrised in terms of the Euclidean 
metric suitably deformed by the Teichmiiller group. This is to say, we can always describe 
(T^jhab) in terms of {T^gi^-^,r]ab), where T]ab is the Euclidean 3d metric and T^.^-^^-^ is the 
deformed torus. The advantage of using the Teichmiiller parameters comes from their easy 
geometrical interpretation in terms of global distortions of a compact manifold, and that 
this scheme can be implemented for each and every possible compact space of a general 
homogeneous (although not necessarily isotropic) metric. 

The specific shape of these distortions is thus going to be parametrised by the Te- 
ichmiiller parameters, that for the 3-torus can be thought of as three-vectors identifying 
the axes of the manifold. Of course we do not need nine parameters to describe the man- 
ifold, as we can always choose the first vector to lie on the, say, x-axis, and the second 
vector to span the xy-plane. In this case the most general deformation can be obtained by 
the following generators 
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acting on the standard Euclidean metric dP = dx^, where ds^ = dt^ — dP. The overall 
volume of the manifold is not fixed, therefore there are six degrees of freedom in the 
Teichmiiller space. The resulting metric is of course 
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The time-dependence of the Teichmiiller parameters can be read off by comparing the 
two metrics (B.2) and (B.4). Of course in simple cases such as the one just worked out, 
it is straightforward (and intuitive) to separate the time dependence from what one calls 
Teichmiiller parameters, in such a way that the metric hij shows explicitly the separation 
in local (time-dependent) and global quantities. 
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As a final example, and because it is relevant for our discussion in the main text, we 
shall specialise further to the usual isotropic FLRW universe, setting equal the three scale 
factors a(t) = b{t) = c{t). If we moreover demand that the angles of the fundamental 
domain are all right, this leaves only diagonal Teichmiiller parameters. The most general 
3d metric for this universe is given by 

dl^ = a{tf {{a\f dx^ + {alf dy^ + (4^ dz^) , (B.5) 

where we have explicitly factored out the scale factor. This compact Bianchi I universe 
has a common expansion rate, but different physical sizes in the three directions of space, 
and therefore will not support an isotropic perfect fluid but rather will have a direction 
dependent equation of state [49, 50]. This means that our parametrisation of the equation 
of state (5.3) is incomplete, and that in general more freedom is allowed when the universe 
is compact. In the most general case there will also be off-diagonal terms describing skew 
toruses, which will carry even more parameters in the attempt to describe the system. 

A final comment on the relevance of the global parameters in the real world. First, as 
it has just been explained, these parameters may impose some restrictions on the specific 
form of the energy-momentum tensor that one can employ; second, as it is well known a 
QFT on a compact manifold is not the same as on a non-compact one, and it is in general 
sensitive to the specific structure of the compact manifold through the boundary conditions 
imposed on the quantum fields. We have given in [6] a concrete example on how the physics 
can be radically different when one quantises on a finite manifold rather than on infinite 
spacetime, where corrections even linear in the size of the fundamental compact domain 
may arise in physical quantities such as the topological susceptibility of the vacuum and 
the value of the chiral condensate. The same phenomenon is expected in Minkowski 4d [1], 
and we have provided several arguments in this work in support of the idea that even an 
expanding universe will exhibit the very same properties. 
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